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The pairwise production of neutral Higgs bosons (h°A°, H°A°) is analyzed in the context of the 
future linear colliders, such as the ILC and CLIC, within the general Two-Higgs-Doublet Model 
(2HDM). The corresponding cross-sections are computed at the one-loop level, including the full 
set of contributions at order 0(al w ) together with the leading 0(a^ w ) terms, in full compliance 
with the current phenomenological bounds and the stringent theoretical constraints inherent to the 
consistency of the model. We uncover regions across the 2HDM parameter space, mainly for low 
tan j3 ~ 1 and moderate - and negative - values of the relevant As parameter, wherein the radiative 
corrections to the Higgs-pair production cross section, <r(e+e- -> A°h /A°H ), can comfortably 
reach |<5(t|/(t ~ 50%. This behavior can be traced back to the enhancement capabilities of the 
trilinear Higgs self-interactions - a trademark feature of the 2HDM, with no counterpart in the 
Minimal Supersymmetric Standard Model (MSSM). The corrections are strongly dependent on the 
actual value of A5 as well as on the Higgs mass spectrum. Interestingly enough, the quantum effects 
are positive for energies around yfs ~ 500 GeV, thereby producing a significant enhancement in the 
expected number of events precisely around the fiducial startup energy of the ILC. The Higgs-pair 
production rates can be substantial, typically a few tens of femtobarn, therefore amounting to a few 
thousand events per lOOfb -1 of integrated luminosity. In contrast, the corrections are negative in 
the highest energy range (viz. y^s ~ 1 TeV and above). We conclude that a precise measurement 
of the 2H final states could carry unambiguous footprints of an extended (non-supersymmetric) 
Higgs sector. Finally, to better assess the scope of these effects, we compare the exclusive pairwise 
production of Higgs bosons with the inclusive gauge boson fusion channels leading to 2H + X final 
states, and also with the exclusive triple Higgs boson production. We find that these multiparticle 
final states can be highly complementary in the overall Higgs bosons search strategy. 

PACS numbers: 12.15.-x, 12.60.Fr, 12.15.Lk 



I. INTRODUCTION 

It was already in the 60 's when the pioneering works of 
Higgs, Kibble, Englert and others suggested the existence 
of a fundamental spinlcss building-block of Nature [l[ , 
whose non- vanishing vacuum expectation value could ex- 
plain the spontaneous breaking of the SU{2)i, <8> U(l)y 
gauge group of the Elcctrowcak (EW) interactions down 
to the U(l) em abelian symmetry group of the Electro- 
magnetism. Indeed, this so-called Higgs mechanism is 
the most fundamental longstanding issue that remains 
experimentally unsettled in Particle Physics. For one 
thing it is the only known strategy capable of building up 
a (perturbative) rcnormalizablc quantum field theoreti- 
cal description of the Elcctrowcak Symmetry Breaking 
(EWSB) phenomenon. It is difficult to overemphasize 
that, to a great extent, the Higgs mechanism embodies 
the backbone of the SM structure and that in its absence 
we would have to cope with an entirely different concep- 
tion of the inner functioning of the SM as a quantum field 
theory (QFT) . 

The central assumption here is the existence of (at 
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least) one elementary spinless field - the so-called Higgs 
boson. Its relevance is particularly evident if we take into 
account that the presence of one or more such fields in 
the structure of the model is indispensable for the unitar- 
ity of the theory. If no Higgs bosons exist below the TeV 
scale, weak interactions would actually become strong at 
that scale and e.g. the W + W~ cross-section would vio- 
late the unitarity of the scattering matrix. Moreover, in 
the absence of Higgs bosons the various particle masses 
could not be generated consistently (i.e. without spoil- 
ing the ultraviolet behavior of the theory at higher orders 
of perturbation theory). In short, if no Higgs bosons are 
there to protect the (presumed) perturbative structure of 
the weak interactions, the latter would enter a perilous 
runaway regime at high energies. It is, thus, essential 
to confirm experimentally the existence of one or more 
Higgs boson particles through their explicit production 
in the colliders. 

Amazingly enough, more than 40 years after the Higgs 
mechanism was naturally blended into the conventionally 
accepted SM landscape of the strong and electroweak in- 
teractions , it remains still thickly curtained and no 
direct evidence has been found yet of its main testable 
prediction: the existence of fundamental spinless parti- 
cles. Notwithstanding the many efforts devoted at LEP 
and at the Tevatron in the last decades, all experimental 
searches for Higgs boson signatures have come up empty- 
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handed, and one has to admit that no conclusive physi- 
cal signature for the existence of these elementary scalar 
fields - not even the single one predicted by the SM - has 
been found for the time being. Thus, we do not know 
whether the conventional Higgs boson exists at all, or if 
there are no Higgs particles of any sort. On the contrary, 
it could be that there are extensions of the SM span- 
ning a richer and elusive spectrum of Higgs bosons which 
have escaped all our dedicated experimental searches. 
Be that as it may, this pressing and highly intriguing 
enigma will hopefully be unveiled soon, specially with 
the advent of the new generation of supercolliders, like 
the brand new LHC, the future International Linear Col- 
lider (ILC) Q and, hopefully, the Compact Linear Col- 
lider (CLIC) too0. 

Therefore, we have to (and we can) be well prepared 
to recognize all kind of hints from Higgs boson physics. 
Let us recall that apart from the single neutral (CP-even) 
Higgs boson H° of the SM, other opportunities arise in 
model building that could serve equally well the afore- 
mentioned purposes. Perhaps the most paradigmatic ex- 
tension of the SM is the Minimal Supersymmetric Stan- 
dard Model (MSSM) Q , whose Higgs sector involves two 
doublets of complex scalar fields. The physical spec- 
trum consists of two charged states, K^, two neutral CP- 
even states h°, H° (with masses conventionally chosen as 
M h o < M H o) and one CP-odd state A . We shall not 
dwell here on the whys and wherefores of supersymme- 
try (SUSY) as a most sought-after realization of physics 
beyond the SM. It suffices to say that it provides a Higgs 
sector which is stable (to all orders in perturbation the- 
ory) under embedding of the low-energy structure into a 
Grand Unified Theory, and in particular it provides the 
most natural link with gravity [5| . 

But, how to test the structure of the Higgs sector? 
In the MSSM case, such structure is highly constrained 
by the underlying supcrsymmctry. A consequence of it, 
which is particularly relevant for our considerations, is 
the fact that the self-interactions of the SUSY Higgs 
bosons turn out to be largely immaterial from the phe- 
nomenological point of view, in the sense that they can- 
not be enhanced (at the tree-level) as compared to the 
ordinary gauge interactions and, therefore, cannot pro- 
vide outstanding signatures of physics beyond the SM. 
The bulk of the enhancing capabilities of the MSSM 
Lagrangian is concentrated, instead, in the rich struc- 
ture of Yukawa-like couplings between Higgs bosons and 
quarks or between Higgs boson and squarks, or even be- 
tween squarks and chargino-neutralinos, as it was shown 
long ago in plentiful phenomenological scenarios involv- 
ing quantum effects on the gauge boson masses @, Q and 
gauge boson and top quark widths @ , as well as in many 
other processes and observables (see e.g. [9l-fl3|. For 
the present state of the art, cf. flil [l5| : and for recent 
comprehensive reviews, see [lH Il7| . for example. 

The two-Higgs S , C/i(2)-doublet structure of the Higgs 
sector in the MSSM is a trademark prediction of SUSY 
invariance. Nonetheless, the very same doublet structure 



can appear in the form of a non-supersymmetric frame- 
work, the so-called general (unconstrained) Two-Higgs- 
Doublet Model(2HDM) , which does also entail a rich phe- 
nomenology [l8[ • Again, the same spectrum of CP-even 
(h°,H°) and CP-odd (A ) scalar fields arise. However, 
we should keep in mind that the most distinctive aspects 
of the general 2HDM Higgs bosons could well be located 
in sectors of the model quite different from the MSSM 
case [lj| . Indeed, even the very couplings involved in the 
structure of the Higgs potential (viz. the so-called trilin- 
ear and quartic Higgs boson self-interactions) could be 
the most significant ones as far as the phenomenologi- 
cal implications are concerned, rather than just the en- 
hanced Yukawa couplings with the heavy quarks. In such 
circumstance, we should be able to detect a very different 
kind of experimental signatures. Assuming, for instance, 
that the LHC will find evidence of a neutral Higgs boson, 
a critical issue will be to discern whether the newcomer 
is compatible with the SM or any of its extensions and, 
in the latter case, to which of these extensions it most 
likely belongs. In order to carry out this "finer" hunting 
of the Higgs boson(s), a TeV-range accessible linear col- 
lider machine (such as the aforesaid ILC/CLIC) will be 
needed. In this work, we will show how useful it could 
be such collider to discriminate neutral supersymmetric 
Higgs bosons from generic (non-SUSY) 2HDM ones. 

The paper is organized as follows. In the next section, 
we remind the reader of some relevant Higgs boson pro- 
duction processes in the linear colliders. In Section Mil 
we describe the structure of the 2HDM and the relevant 
phenomenological restrictions, including the properties 
of the trilinear Higgs self-interactions and their interplay 
with the notion of unitarity and vacuum stability. In 
Section ITVl we develop in detail the renormalization pro- 
gram for the Higgs sector of the general 2HDM. The the- 
oretical setup for the one-loop computation of the Higgs 
boson cross-sections is elaborated in Section [V] leaving 
Section lVTl to present a comprehensive numerical analysis 
of our results. Finally, Section I VIII is devoted to a thor- 
ough general discussion and to present the conclusions of 
our work. 



II. HIGGS BOSON PRODUCTION AT LINEAR 
COLLIDERS 

The leading Higgs boson production mechanisms at the 
LHC have been studied thoroughly for the last twenty 
years, they are well under control both in the SM [2(| 
and in the MSSM [U [22[ and there are also some studies 
in the 2HDM [23j . Hopefully, they may soon help reveal- 
ing some clues on Higgs boson physics at the LHC - 
for a review see e.g. [171 ] and references therein. How- 
ever, it is not obvious that it will be possible to easily 
disentangle the nature of the potentially produced Higgs 
boson(s). The next generation of TeV-class linear collid- 
ers (based on both e + e~ and 77 collisions), such as the 
ILC and the CLIC projects [1,0], will be of paramount 
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importance in finally settling the experimental basis of 
the "Higgs issue" as the most fundamental theoretical 
construct of the SM of electroweak interactions. Thanks 
to its extremely clean environment (in contrast to the 
large QCD background inherent to a hadronic machine 
such as the LHC), a linear collider (linac for short) should 
allow for a precise measurement of the Higgs boson pa- 
rameters, such as: i) the Higgs boson mass (or masses, if 
more than one); ii) the couplings of the Higgs bosons to 
quarks, leptons and gauge bosons; iii) the Higgs boson 
self-couplings mentioned above, i.e. the trilinear (3H) 
and quartic (4H) Higgs boson self-interactions. At the 
end of the day, a linac should allow us to dig deeper than 
ever into the structure of the EWSB and, hopefully, even 
to reconstruct the Higgs potential itself. 

A detailed road map of predictions for Higgs-boson ob- 
servables at the linear colliders is called for, with a special 
emphasis on those signatures which may be characteris- 
tic of the different extensions of the SM. For example, 
as indicated above, triple Higgs boson self-interactions 
(3H) may play a cardinal role in this endeavor because, 
in favorable circumstances, they could easily distinguish 
between the MSSM and the general 2HDM Higgs sec- 
tors. Such 3H couplings can mediate a plethora of in- 
teresting processes. The key point here is the potentially 
large enhancements that the 3H couplings may accommo- 
date in the general 2HDM case, in contrast to the super- 
symmetric extensions. Actually, even for the SUSY case 
there is a large number of works attempting to extract 
vestiges of non-standard dynamics in these couplings, 
mainly through the radiative corrections that they can 
undergo. For instance, the 3H-couplings have been in- 
vestigated in (23 - |26| . and in some cases considered for 
possible phenomenological applications in TeV-class lin- 
ear colliders through the double-Higgs strahlung process 
e + e~ — > KHZ or the double-Higgs WW-fusion mecha- 
nism c + c~ — >• H + H~^ c 77 c . These processes, which include 
vertices like ZZH, WWH, ZZHH, WWHH and HHH, are 
possible both in the SM and its extensions, such as the 
MSSM and the general 2HDM. Unfortunately, the cross- 
section turns out to be rather small both in the SM and 
in the MSSM, being of order 10 -3 pb at most, i.e. equal 
or less than 1 fb |25[ ■ Even worse is the situation regard- 
ing the triple Higgs boson production in the MSSM, in 
which - except in the case of some particular resonant 
configuration - the typical cross-sections just border the 
line of ~ 0.01 fb or less [25j . In the latter reference, for 
instance, it has been shown that if the double and triple 
Higgs production cross-sections would yield sufficiently 
high signal rates, the system of couplings could in prin- 
ciple be solved for all trilinear Higgs self-couplings up 
to discrete ambiguities using only these processes. But 
this is perhaps a bit too optimistic since in practice these 
cross-sections are manifestly too small to be all measur- 
able in a comfortable way. 

In stark contrast to this meager panorama within the 
SM and the MSSM, 3H couplings have been tested for 
tree-level processes in the context of the unconstrained 



2HDM. For example, in Ref. [27} the tree-level produc- 
tion of triple Higgs boson final states was considered in 
the ILC. There are three classes of processes of this kind 
compatible with CP-conservation, namely 

1) c + c~ -> H + H~h , 2) e + c~ -> hhA° , 

3) c+c- -> h°H°A , (h = h°, H°, A ) (1) 

where, in class 2), we assume that the two neutral Higgs 
bosons h must be the same, i.e. the allowed final states 
are (hhA°) = (h°h°A°), (H°H°A ) and (A A A ). The 
cross-sections in the 2HDM were shown to reach up to 
0(0.1) pb [l?}, i.e. several orders of magnitude over the 
corresponding MSSM predictions (28|. Besides the exclu- 
sive production of three Higgs bosons in the final state, 
also the inclusive pairwise production of Higgs bosons 
may be critically sensitive to the 3H self-interactions. 
Sizable production rates, again in the range of 0.1 — lpb 
have recently been reported in Ref. [29[ , whose focus was 
on the inclusive Higgs boson-pair production at order 
®( a tw) through the mechanism of gauge-boson fusion 

e + c~ -> V*V* -> hh + X 

{V^W ± ,Z, /i = h°,H°,A°,H ± ). (2) 

Similar effects have also been recently computed on 2H- 
strahlung processes of the guise e + e~ — > Z° h h [30[ • In a 
complementary way, 3H couplings can also be p robed in 
loop-induced processes. For instance, Ref. [3l[ presents 
a computation of the single Higgs boson production rate 
through the scattering process of i) two real photons, us- 
ing the 77 mode of a linear collider, i.e. 77 —¥ h; and ii) 
the more traditional mechanism of virtual photon fusion 
in c + e~ colliders, e + c~ — > 7*7* —th + X. In cither case, 
the obtained cross sections within the general 2HDM (up 
to lpb for the first mechanism, and 0.01 pb for the lat- 
ter) are 1 to 2 orders of magnitude larger than the ex- 
pected SM yields and also well above the MSSM results 
(see [31I ] and references therein). Promising signatures 
have also been reported within the general 2HDM us- 
ing loop-induced production of two neutral Higgs bosons 
through real 77 collisions [32[, although in this case the 
cross-sections are smaller than in the primary mechanism 
77 — > h mentioned before. 

A crucial observation concerning our aim here is the 
following: all the processes described above are directly 
sensitive to the 3H self-couplings already at the leading 
order. In this paper, we continue exploiting the prop- 
erties of these couplings in the general 2HDM, but we 
concentrate now on their impact at the level of indirect 
effects through radiative corrections. The latter may sig- 
nificantly affect processes that are kinematically more 
favored (e.g. because they have a smaller number of 
particles in the final state), but which are nevertheless 
totally insensitive to the trilinear Higgs boson couplings 
at the lowest order. Specifically, we wish to concentrate 
on identifying the largest quantum effects that the 3H 
couplings can stamp on the cross-sections for two-body 
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neutral Higgs boson final states: 



e + e 



2h 



(2^h A ;H A ^ 



(3) 



Surprisingly enough, very little attention has been paid 
to these basic production processes and to the calcula- 
tion of the corresponding radiative corrections within the 
2HDM. In contrast, a lot of work has been invested on 
them from the point of view of the MSSM H [H, |H- 
HH - see also [39[ and references therein 1 . To the best 
of our knowledge, only Refs. j4lTl43j have addressed this 
topic in the general 2HDM, although they are restricted 
to the production of charged Higgs pairs H + H~. In the 
present paper, we shall be concerned exclusively on com- 
puting the quantum effects involved on the neutral Higgs 
boson channels ([3]) at order 0(a^ w ) (and eventually also 
the leading 0(a^, w ) terms). As we will see, very large 
quantum effects (of order 50%) may arise on the cross- 
sections of the two-body Higgs boson channels (j3|) as a 
result of the enhancement capabilities of the triple (and, 
to a lesser extent, also the quartic) Higgs boson self- 
interactions. These effects are completely unmatched 
within the MSSM and should therefore be highly charac- 
teristic of non-supcrsymmetric Higgs boson physics. 



III. THE TWO-HIGGS-DOUBLET MODEL: 
GENERAL SETTINGS AND RELEVANT 
RESTRICTIONS 

The Two-Higgs-Doublet Model (2HDM) is defined 
upon the canonical extension of the SM Higgs sector with 
a second S[/(2)l doublet with weak hypcrchargc Y=l, 
so that it contains 4 complex scalar fields arranged as 
follows: 



$i = 



<f>+ 



(Y = +l), *2=(J 2 oJ (Y = +l). 

(4) 

With the help of these two weak-isospin doublets, we can 
write down the most general structure of the Higgs po- 
tential fulfilling the conditions of CP-conservation, gauge 
invariance and renormalizability. Moreover, a discrete 
Z 2 symmetry $, (-!)*$, (i = 1,2) - which will 
be exact up to soft-breaking terms of dimension 2 - is 
usually imposed as a sufficient condition to guarantee 
the proper suppression of the Flavor-Changing Neutral- 
Current (FCNC) effects that would otherwise arise within 
the quark Yukawa sector 2 (we shall return to this 



point later) 3 . All in all, one arrives at the following ex- 
pression for the tree-level potential: 



2s 2 
2 



+ A 3 

+A 4 
+A 5 
+A 6 



($t*i)($5$ 2 ) - ($t* 2 )(*£$0 



-i 2 



Ke($ t 1 $ 2 ) - 
3 m (${$ 2 ) 



vi v 2 



(5) 



where Aj (i = 1, . . . 6) are dimensionless real parameters 
and Vi (i = 1,2) stand for the non- vanishing VEV's that 
the neutral component of each doublet acquires, normal- 
ized as follows: ($") = v t /V2. 

The complex degrees of freedom encoded within each 
Higgs doublet in Eq.Q are conveniently split into real 
ones in the following way: 



$ 2 = 



2 - 



02 + 
V2 



(0) 



Upon diagonalization of the Higgs potential ([5]) we may 
obtain the physical Higgs eigenstates in terms of the 
gauge (weak-eigenstate) basis: 



H° 
h° 

G° 
A" 

G ± 
H± 



cos a sin a 

— sin a cos a 

cos /3 sin j3 

— sin f3 cos j3 

cos (3 sin (3 
- sin (3 cos (3 



(7) 



The parameter tan/3 is given by the ratio of the two 
VEV's giving masses to the up- and down-like quarks: 



tan/3 



(*?> 



V2 
''1 



(8) 



Incidentally, let us note that this parameter could be 



ideally measured in e + e colliders [36( since the charged 



Higgs event rates with mixed decays, e + e~ — > H H~ — > 
tbrP r , thv T f , do not involve the mixing angle a. The 
physical content of the 2HDM embodies one pair of CP- 
even Higgs bosons (/, H°); a single CP-odd Higgs boson 



1 For the analysis of the tree-level double Higgs production pro- 
cesses in the MSSM, see e.g. the exhaustive overview |40H . For 
a detailed 'anatomy' of the MSSM Higgs sector and an updated 
account of its phenomenological consequences, see e.g. 11711 ■ 

2 For alternative strategies of building up realizations of the 2HDM 
with no explicit Z2 symmetry, see e.g. 14511 and references therein. 



3 This symmetry is automatic in the MSSM case, although it 
is again violated after introducing the dimension 2 soft-SUSY 
breaking terms. These are essential for the EWSB, which other- 
wise would not occur, see Eq. H15I I below. 
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(A ); and a couple of charged Higgs bosons (H*). The 
free parameters of the general 2HDM are usually chosen 
to be: the masses of the physical Higgs particles; the 
aforementioned parameter tan /3; the mixing angle a be- 
tween the two CP-even states; and finally the coupling A5, 
which cannot be absorbed in the previous quantities and 
becomes tied to the structure of the Higgs self-couplings. 
To summarize, the vector of free inputs reads 



(M^o , Mho , Mjsfi , Mh± , a, tan /3, A5) . 



(9) 



Therefore, we are left with 7 free parameters, which in- 
deed correspond to the original 6 couplings A^ and the 
two VEV's vi , V2 - the latter being submitted to the con- 
straint v 2 = v\ + v\ — 4 M^y/g 2 , which is valid order by 
order in perturbation theory, where Myy is the W mass 
and g is the SUl{2) gauge coupling constant. The di- 
mension 2 term that softly breaks the Zi symmetry can 



be written as 



'12 



$f $2 + h.c, with 

Gp 1 tan/3 
2 sin/3cos/3 = F 



2V2 1 + tan 2 P 



the second equality being valid at the tree-level. Given 
tan /3 one may trade the parameter A5 for mi 2 through 
this formula, if desired. 

Since we are goin g to adopt the on-shell renormaliza- 
tion scheme [461-4 8 1 ] for the one-loop calculation of the 
cross-sections (cf. section |PV]) . it is convenient to intro- 
duce the electromagnetic fine structure constant a em = 
e 2 /47r, which is one of the fundamental inputs in this 
scheme (together with the physical Z° mass, Mz, and the 
Higgs and fermion particle masses). The electron charge 
— e is connected to the original SUl(2) gauge coupling 
and weak mixing angle 6w through the well-known rela- 
tion e = g sin 9w , which is also preserved order by order 
in perturbation theory. Then the W ± mass can be re- 
lated to a em and the Fermi constant Gf m the standard 
way in the on-shell scheme, namely 



Gf 

7f 



24 



(1 - Ar) . 



(11) 



where the parameter Ar [48[ vanishes at the tree-level, 

but is affected by the radiative corrections to /x-decay 
both from standard as well as from new physics, for in- 
stance from the MSSM @, and also from the 2HDM 
(see below). Notice that, in the above formula, My/ 
also enters implicitly through the relation sin 2 9\y = 
1 — Alyy/Mg , which is valid order by order in the on- 
shell scheme. Since Gf can be accurately determined 
from the /z-decay, and Mz has been measured with high 
precision at LEP, it is natural to take them both as ex- 
perimental inputs. Then, with the help of the non-trivial 
relation (fTT]) , the W ± mass can be accurately predicted in 
a modified on-shell scheme where Gf, rather than Myy, 
enters as a physical input. 

We shall not dwell here on the detailed structure of the 
Yukawa couplings of the Higgs boson to fermions (see 



[T§ | for an exhaustive treatment), since their contribu- 
tion to the processes e + e~ — > h Q A°/H A under study 
at the one-loop level is largely subdominant in the main 
regions of parameter space. However, it is precisely the 
coupling pattern to fermions that motivates the distinc- 
tion between the different types of 2HDM's, so let us 
briefly describe them. Assuming natural flavor conserva- 
tion [44| , we must impose that at most one Higgs doublet 
can couple to any particular fermion type. As a result, 
the coupling of the 2HDM Higgs bosons with fermions 
(say the quarks) can be implemented in essentially two 
manners which insure the absence of potentially danger- 
ous (tree-level) FCNC's, to wit (cf. Table |Inj! 4 : i) type-I 
2HDM, wherein the Higgs doublet ($2) couples to all of 
quarks, whereas the other one ($1) does not couple to 
them at all; or ii) type-II 2HDM, in which the doublet 
<3?i(resp. $2) couples only to down-like (resp. up-like) 
right-handed quarks. In the latter case, an additional 
discrete symmetry involving the chiral components of the 
fermion sector must be imposed in order to banish the 



A5 j (10) tree-level FCNC processes, e.g. D\ 



U'r. 



for the down and up right-handed quarks in the three fla- 
vor families (i = 1, 2, 3). 

In turn, the Higgs self-couplings A; in the Higgs poten- 
tial can be rewritten in terms of tan /3 and the physical 
parameters of the on-shell scheme such as the masses and 
the electromagnetic fine structure constant a em . At the 
tree-level, we have 



Ai = 



A 5 (1 — tan 2 



2M 2 v s 2 w cos 2 /3 



[M% cos 2 a + Mia sin 2 a 
-(Mfjo - M 2 ) sin 2a cot/3 



A 5 (1 - 1/tan 2 P) 



4 o A/f2 ,2 „;„2 

f2 „„„2 „, 1 7i,r2 _;_2 



2 M 2 , s- w s: 



x [M 2 cos 2 a + Mfjo sin 2 a 
--(Mjjo - Ml a ) sin 2a tan/3 



A3 = - — 



Ao 
4 

2 OLpir 



M 2 s 2 
lvl w b w 



a em 7r sin 2a 
2M 2 v s 2 w sin 2p 

■ M 2 H± , 



(M 2 H0 



M 2 ). 



A fi = 



2a e 



M 



M 2 A a, 



(12) 

cos 6w From these La- 
grangian couplings in the Higgs potential, we can de- 



whe 



re sw 



w °w 
sin 6w, cw 



4 Strictly speaking, in the absence of CP-violation there are four 
manners to have natural flavor conservation, giving rise to four 
different kinds of allowed 2HDM's that do not lead to tree-level 
FCNC's, see [49l |. although we will restrict ourselves here to the 
two canonical ones llSH . 
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type I 


type II 


11 Li 


cos ol 1 sin (3 


cos a 1 sin P 


h°hh 
11 uu 


CUb (_i / bill 


— bill (JC 1 LUb fj 


H u tt 


sin a/ sin p 


sin a/ sin f) 


E%b 


sin a/ sin p 


cos a / cos P 


A u tt 


cot /? 


cot/3 


A°66 


-cot/9 


tan /3 



TABLE I: Neutral Higgs boson couplings to fermions in type-I 
and type-II 2HDM, using third family notation. For exam- 
ple, the SM Higgs boson coupling to top and bottom quarks, 
—gmf/2Mw (/ = t,b), multiplied by the corresponding fac- 
tor in the table provides the 2HDM couplings of the h° boson 
to these quarks. Worth noticing are the characteristic en- 
hancement factors arising at large (and low) values of tan p. 



rive the "physical couplings" , namely those affecting the 
physical Higgs bosons in the mass-eigenstate basis. We 
call these couplings the triple (3H) and quartic (4H) 
Higgs couplings. Their behavior and enhancement ca- 
pabilities are at the very core of our discussion. Indeed, 
from our analysis it will become clear that they furnish 
the dominant source of quantum corrections to the Higgs- 
pair production processes © within the general 2HDM. 

The physical 3H and 4H couplings are not explicitly 
present in the 2HDM potential (0. They are derived 
from it after spontaneous breaking of the EW symme- 
try and corresponding diagonalization of the Higgs boson 
mass matrix using the rotation angles a and /3 in ([7]) . In 
the particular case of the SM, the trilinear and quartic 
Higgs couplings have fixed values, which depend uniquely 
on the actual mass of the Higgs boson. In the MSSM, 
however, and due to the SUSY invariance, the Higgs bo- 
son self-couplings are of pure gauge nature, as we shall 
revise briefly below. This is in fact the primary reason for 
the tiny production rates obtained for the triple Higgs bo- 
son processes ([T]) within the framework of the MSSM [25| . 
In contrast, the general 2HDM accommodates trilinear 
Higgs couplings with great potential enhancement. The 
full list is displayed in Table [H] 

As can be seen, the couplings in this table depend 
on the 7 free parameters ([5]). In the particular case 
where A 5 = A 6 , this table reduces to Table 1 of Ref. [27| ■ 
The equality of these couplings takes place automatically 
e.g. in the MSSM, where in addition other simplifica- 
tions occur, as we discuss below. Let us note, for exam- 
ple, that in the limit a = f3 — tt/2 the h°h°h -trilinear 
coupling in Table [TT] reduces exactly to the SM form 
^hhh = ~3igM%/(2M w ) - where we denote by H 
the SM Higgs boson. This situation would correspond to 
the so-called decoupling limit in the MSSM [f3(| since it 
is correlated with Mao — > oo, although there is no such 
correlation in the general 2HDM. 

Let us recall that the MSSM Higgs sector is a type-II 
one of a very restricted sort: it is enforced to be SUSY 
invariant. This is a very demanding requirement as, for 
example, it requires that the superpotcntial has to be a 



holomorphic function of the chiral supcrficlds [H, [HI). 
This implies that we cannot construct the MSSM with 
two Y = +1 Higgs supcrficlds. Then, in order to be able 
to generate masses for both the top and bottom quarks 
through EWSB, the $2 doublet can be kept as it is with 
Y = +1, although we must replace $1 with the conjugate 
(Y = — 1) SU L (2) doublet 

) Se *;=(5_) (Y=-D, (13) 

where e = i<72, with <ji the second Pauli matrix. Thus, 

for the first doublet the correspondence with the MSSM 
case reads $1 = — effj\ and the second doublet (the one 
with no change) is just relabeled i?2- Moreover, the tree- 
level relationships imposed by Supersymmetry between 
the Ai couplings of the potential ([5]) are the following: 

Ai = A2 , 



_ 2-na em 

A 4 — LA\ ^ > 

C W 

■\ \ r> \ 2tt ty ern f . 

A 5 = A 6 = 2Ai - — — — . (14) 

s w c w 

Substituting (fl~3|) and (|14| in ((5]) one obtains the usual 

MSSM potential, which in practice must be supple- 
mented with soft SUSY-breaking scalar mass terms 
mf Hf , including a bilinear mixing term m\ 2 H± H2 for 
the two Higgs doublets. This term is the analog of (fT0|) 
for the 2HDM, but in the soft-SUSY breaking context is 
arbitrary. The result can be cast as follows: 

V(H u H 2 ) = (\rf+ml) iH^ + Hrf+mi) \H 2 \ 2 

zs w c w s w 
-m? a (e«flj H{ + /i.e.) , (15) 

where \i is the higgsino mass term in the superpotcn- 
tial Q . The obtained potential is one where all quartic 
couplings become proportional to a em , which is tanta- 
mount to say that in the SUSY case these self-couplings 
arc purely gauge. Thus, after EWSB also the trilinear 
self-couplings will be purely gauge. As warned before, 
this is the primary reason for their phenomcnological in- 
conspicuousncss. Under the supcrsymmctric constraints, 
the entries of TableHQboil down to the MSSM form listed 
e.g. in [l8|. Moreover, the five constraints (fTl|) reduce 
the number of free parameters from 7 to 2 in the su- 
persymmetric context, typically chosen to be tan/3 and 
M40. For the general 2HDM case, however, we shall stick 
to the form presented in Table HIl where the 7 free inputs 
are chosen as in ©. 
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h°h°h° 


~2m w ^20s w i M H'> ( 2 co < a + P) + sin 2 « sin(/? - a)) 
cos(a + /J)cos 2 (/3 a) 4X5 M \ S H 


h°h°H° 


-2M™^L + M^ ) sin2a - (3sin 2a - sin 2/3) 


h°H°H° 


™T£wl w [(A4o+2A4 )sin2a-(3sin2a + sin2/3)^ 2*L*&] 


H°H°H° 


[M| (2sin(a + /?) - cos(/J - a) sin 2a) 
-sin(a + /J)sin 2 (/3-a) S ^ i^^l 


h°A°A 


ie [cos( Q +/3) / OJ , f 2 «5< s2 w) „■ iq \ ( M 2 oj\f2 \] 


n A Cj 


2M w s w \ M A» M h°) COS (P ") 


h°G°G° 


- 2A / e s Af 2 o sin^-a) 


H°A (, A° 


ie r s in( Q + /3) (c, M 2 4A 5 «w 4 \ po „fa „.\ /« f 2 r "\fMl 
2M W S W [ sin2/3 ^"ffU e 2 J "MP aj^VJ^o ZIV1 a q)^ 


ttO \O r <0 
II \ ( j 


ie ( \ j 2 u2 \„;„/o „ \ 

2M WSW - l M A° A%o) sm (P a) 


H°G°G° 


-2Mw S „- M ho cos(/3-a) 


h°H+H- 


- wfe [^W 1 ( 2M *° " ^f^) - (M 2 o - 2M£_) sm(/3 - a)] 


H°H+ET 


-ni£nF ( 2M -° " ^#^) - cos(/3 - a) (A4o - 2M^_ )] 



TABLE II: Trilinear Higgs boson self- interactions (iXm) in the Feynman gauge within the 2HDM. Here G° is the neutral 
Goldstone boson. These vertices are involved in the radiative corrections to the A h° Z°(A° H° Z°) bare couplings. 



In another vein, it is of paramount importance when 
studying possible sources of unconventional physics to 
make sure that the SM behavior is sufficiently well re- 
produced up to the energies explored so far. Such a re- 
quirement translates into a number of constraints over 
the parameter space of the given model. In particular, 
these constraints severely limit the a priori enhancement 
possibilities of the Higgs boson self-interactions in the 
2HDM. 

• To begin with, we obviously need to keep track of 
the exclusion bounds from direct searches at LEP 
and the Tevatron. These amount to Mi > 114GcV 
for a SM-like Higgs boson - although the bound is 



weakened down to 92.8 GeV for models with more 
than one Higgs doublet [Ht In triC 2HDM all the 
mass bounds can be easily satisfied upon a proper 
choice of the different Higgs masses - which, unlike 
the SUSY case, are not related with each other. 

• Second, an important requirement to be enforced 
is related to the (approximate) SUl+r(2) custo- 
dial symmetry satisfied by models with an arbi- 
trary number of Higgs doublets (52|. I n practice, 
this restriction is implemented in terms of the pa- 
rameter p, which defines the ratio of the ncutral- 
to-charged current Fermi constants. In general it 
takes the form p = p$ + Sp, where po is the tree- 
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level value. In any model containing an arbitrary 
number of doublets (in particular the 2HDM), we 
have po = M^/MgC^ = 1, and then Sp represents 
the deviations from 1 induced by pure quantum 
corrections. From the known SM contribution and 
the experimental constraints [5ll | we must demand 
that the additional quantum effects coming from 
2HDM dynamics satisfy the approximate condition 
\$P2HDm\ ^$ 10~ 3 . It is thus important to stay in 
a region of parameter space where this bound is 
respected. In our calculation we include the domi- 
nant part of these corrections, which involves one- 
loop contributions to Sp mediated by the Higgs- 
boson and yields [53| 



dp 



G, 



2HDM = 



8^ 



M 2 H± 



1 



M 2 
M A o 



+ cos 2 (J3 - a)M l 



h° 



M 2 A0 



m^ 1 



M 2 h o 



M 2 



In- 



M 



2 1 



Ml 



+ sin 2 (fi -a) M 2 m 



M 2 

1VI A0 



M 2 



M 2 



In 



M 2 H± - M 2 H , 



In 



M 



M 2 
M H° 



(16) 



From this expression, it is clear that arbitrary mass 
splittings between the Higgs bosons could easily 
overshoot the limits on 5p. However, we note that 
if M A o — > M H ± then <5p2HDM — > 0, and hence if the 
mass splitting between M A o and M H ± is not signi- 
ficant (5/02HDM can be kept under control 5 . Let us 
also recall in passing that the Sp correction trans- 
lates into a contribution to the parameter Ar of 
(ITT1) given by — {c 2 v /s 2 v )8p. Therefore, the tight 
bounds on Ar, together with the direct bounds on 
the ratio of charged and neutral weak interactions, 
restrict this contribution as indicated above. 

Also remarkable are the restrictions over the 
charged Higgs masses coming from FCNC radia- 
tive B-meson decays, whose branching ratio B(b — > 
sj) ~ 3 x 10 -4 |51| is measured with sufficient 
precision to be sensitive to new physics. Current 
lower bounds render Af H ± > 295 GeV for tan /3 > 1 
[55l |5(| . It must be recalled that these bounds ap- 
ply for both Type-I and Type-II models in the re- 
gion of tan/3 < 1. In contrast, they are only rel- 
evant for Type-I in the large tan /3 regime, since 
for them the charged Higgs couplings to fermions 
are proportional to cot /3 and hence the loop con- 
tributions are highly suppressed at large tan j3 - cf. 



Ref.[57| for detailed analytical expressions of the 
contributions to B(b — > 57) from the 2HDM. Be- 
sides, data from B — > Ivi may supply additional re- 
strictions on the charged Higgs mass at large tan f3 
for the type-II 2HDM (58|. For a very recent up- 
dated analysis of the different i?-physics constraints 
on the 2HDM parameter space, see Ref. (45| . 

• Further constraints apply to tan /3 coming from the 
following two sources: i) The Z° — > bb and B — B 
mixing processes strongly disfavor tan/3 below 1 
56| ; ii) The requirement that the Higgs couplings 

to heavy quarks remain perturbative translates 
into an (approximate) allowed range of O(0. 1) < 
tan/3 < 60 j49|. 

• Besides the available experimental data, an addi- 
tional set of requirements ensues from the theoreti- 
cal consistency of the model. In particular, we shall 
introduce the following set of conditions ensuring 
the stability of the vacuum [59Tj6l| : 



Ai + A 3 > ; 
2 v / (Ai+A 3 )(A 2 
+min[0, A 5 — A 4 , A 6 



A 2 + A 3 > 
' A 3 ) + 2A 3 4 
A 4 ] > 0. 



A 4 



(17) 



The unitarity bounds deserve a separate discussion. In- 
deed, they turn out to be a fundamental ingredient of our 
computation. As we have already been mentioning, the 
trademark behavior that we aim to explore shall depend 
critically on the enhancement capabilities of the Higgs- 
boson self-interactions which, in turn, will be sharply 
constrained by unitarity. The basic idea here is that, 
within perturbative QFT, the scattering amplitudes are 
"asymptotically flat" , meaning that they cannot grow in- 
definitely with the energy. This is tantamount to say that 
the unitarity of the S-matrix must be guaranteed at the 
perturbative level. In a pioneering work by Lee, Quigg 
and Thacker (LQT) (62|, the above arguments were first 
applied to the SM. As a very relevant outcome, an upper 
bound for the Higgs boson was obtained: 



Ml < 8nv 2 = 



8 71 



^ = M LQT^( L2TeV ) • ( 18 ) 



The basic idea beneath the LQT argument is to com- 
pute the scattering amplitudes in a variety of processes 
(namely scalar-scalar, scalar-vector and vector-vector 
scattering; most particularly, the longitudinal compo- 
nents of the vector bosons) and demand all of them to be 
in accordance with the generic tree-level unitarity condi- 
tion. The latter is a pure quantum-mechanical require- 
ment, which can be expressed in terms of the s-wave scat- 
tering lenght in the high-energy limit: 



M < 1/2. 



(19) 



5 For more refined analyses on Sp constraints, see e.g. [54| . 



Taking into account the equivalence between the longitu- 
dinal components of the gauge and the Goldstone bosons 



9 



(which is exact to order 0(yj Myy/ s), s being the center- 
of-mass energy squared), we can convince ourselves that 
only the scalar-scalar processes will be relevant to this 
concern. Moreover, it can also be shown that the bulk of 
the contribution is brought about by the quartic scalar 
self-couplings - the triple ones are suppressed as 1/s. 
Thereby the quartic Higgs self-couplings (4H) turn out 
to be the most strongly constrained by the bounds of 
Eq. ™. 



h 
Pi 



J_ {2(A3+A4) _^ } 



(20) 



IV. RENORMALIZATION OF THE 2HDM 
HIGGS SECTOR 



Several authors have exported the above ideas to the 
2HDM I63l-l76ll and also to its complex (CP- violating) ex- 
tensions [77j . The strategy here is to compute the S- 
matrix elements Sij in a set of 2 — > 2 scattering pro- 
cesses involving Higgs and Goldstone bosons. In this 
case there are more scalar-scalar channels than in the 
SM, which are coupled among themselves. At this point, 
we could perform a unitary transformation U relating 
the S-matrix expressed in the original weak-interaction 
cigenstatc basis (S w ) with the S-matrix written in the 
physical mass-eigenstate basis (S m ). Since, however, the 
former is much simpler than the latter and both carry ex- 
actly the same information, we may generalize the origi- 
nal LQT procedure |62| by focusing on S w and requiring 
that its eigenvalues satisfy the tree-level unitarity condi- 
tion |oj| < 1 (Vi). In practice, inspired by the scattering 
length analysis (|T9|) . we require \on\ < 1/2 (Vi). This 
leads to a set of conditions over several linear combina- 
tions of the quartic couplings in the original Higgs po- 
tential 68 76], as can be seen from the explicit form of 
the various eigenvalues on at the tree-level: 



a± 



7^{3(A 1+ A 2 + 2A 3 ) 



±(j9(Ai - A 2 ) 2 + (4A 3 + A 4 + y + ^) 2 } 



b± = ^{Ai + A 2 + 2A 3 



- a 5 + \ 6 y j 



±A /(A 1 -A 2)2 + ^ A4 + 
c± = d± = — {Ai+A 2 + 2A 3 

167T I 



±\/(Ai-A 2 ) 2 



(a 5 - x 6 y 



A5 Ag 
~2 2 

A.5 Ag 

Y ~ ~2 



6ni' 

f+ = -^{2A 3 -A 4 + 5 

167T I 

/- = T7^f 2A 3 + A 4 

107T I 



}■ 
}■ 
}■ 
}■ 
}■ 



In this Section, we discuss in detail the rcnormaliza- 
tion of the Higgs sector. The renormalization of the SM 
fields and parameters is performed in the conventional 
on-shell scheme in the Feynman gauge, see e.g. Ref. 
(46M481 [zi, [zi|. At present, highly automatized proce- 
dures are available for loop calculations, especially in the 
MSSM, see e.g. [l6[ and [80l- 83] . However, in our case 
the calculation of the cross-sections for the processes Q 
is performed within the general (non-supcrsymmetric) 
2HDM and we must deal with the renormalization of the 
Higgs sector in this class of generic models. To this end, 
we attach a multiplicative wave- function (WF) renormal- 
ization constant to each of the SUl{2) Higgs doublets in 
the 2HDM. 



$0 



-1/2 / $+ 



<*> 2 + 



I $0 



(21) 

The renormalizcd fields arc those on the r.h.s. of these 
expressions. At one-loop we decompose Z§ 1 2 = 1 + 
5Zq, 1 2 + 0(a ew ). These WF renormalization constants 
in the weak-eigenstate basis can be used to construct the 
WF renormalization constants Z^h- = 1 + 8Z} ii h i m the 
mass-eigenstate basis by means of the set of linear trans- 
formations f7|). We are interested only in the production 
of neutral Higgs bosons in this paper, and therefore we 
shall only quote the relations referring to them. The cor- 
responding 8 Zhihj in the neutral Higgs sector (using the 
notation 8 Zh t hi = $ %hi) are determined as follows: 

8 Z\fi = sin 2 a 8 Z$ 1 + cos 2 a 8 Z$ 2 

8 2 H o = cos 2 a 8 Z^ 1 + sin 2 a 8 Z$ 2 

8 Z A o = sin 2 /? 8 Z<s> x + cos 2 (3 8 Z^ 2 

8 Z G o = cos 2 (3 8 Z<t> 1 + sin 2 /35 Z$ 2 



8 Zho-fio = sin a cos a (8 Z$, 2 — 8 Z^, 1 ) 
(5 Zap g o = sin j3 cos j3 {8 Z^ 2 — 8 Z<$> 1 ) , 



(22) 



where in the last two equations we have also included 
the mixing of the two CP-even Higgs bosons and that of 
the CP-odd state with the neutral Goldstone boson in 
the Feynman gauge. The different counterterms 8 Z^h 
must be anchored by specifying a set of subtraction con- 
ditions on a finite number of 2-point functions. To that 
purpose, we need to specify the renormalizcd self-energies 
corresponding to the gauge bosons and each of the phys- 
ical Higgs-boson fields in the model. First of all let us 
define the relation between the polarization tensor of a 
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gauge boson and the transverse self-energies. We define 
it as follows, 



(q 2 ) 



(23) 



where the dots indicate that we neglect the longitudinal 
part (proportional to q^q v ) because our calculations for 
a linear collider involve extremely light external fcrmions 
(electrons and positrons) and as a result the longitudinal 
contributions are suppressed as m 2 /M v . The same ap- 
plies to the q^q" terms in the transverse part, of course. 
We identify +iT, fll/ (q) with the Feynman diagram asso- 
ciated to the structure of the polarization tensor. As a 
result, the unrenormalized self-energy part is defined in 
such a way that — iT,vv(q 2 ) (notice the minus sign) is 
equated to the blob of external lines V and V. For sim- 
plicity, we define Ey(g 2 ) = Eyy(g 2 ) and similarly for 
other fields. Thus e.g. the free propagator for the gauge 
boson V becomes "dressed" by the quantum effects as 
follows: 



M. 



(0)2 



q 2 - M v 



(0)2 



Zv(q 2 



-igp.u z v 



q 2 -M v + E y (g 2 



(24) 



My' being the bare mass of V and where in the sec- 
ond equality we have introduced the renormalizcd form 
of the propagator. Similarly, for self-energies involving 
two scalars (cf. Fig.[T]) we identify +iY l( f )i< f, j (q 2 ) with the 
Feynman blob with external legs <pi and <f>j. For instance, 
the dressed propagator of the scalar <j> with bare mass 

(0) • u 

m\ is given by 



i Za 



q 2 - ?4 0)2 + E (g 2 ) 



q 2 - m 2 + E (g 2 ) 



(25) 



The last equality refers to the corresponding renormal- 
ized scalar field propagator, being the renormal- 
ized mass and E^, the rcnormalized self-energy. Without 
wanting to appear too pedagogical, we remind the reader 
that the renormalization transformation at one-loop is 
performed by canonically introducing renormalization 
constants relating every bare field 4>i to the renormal- 



izcd ones as follows 



,(0) 



,1/2 



= (1 + \&Zi)<l>i - 
which in this case takes on the particular form (|21l) 
and at the same time decomposing the bare masses and 
couplings as the sum of the rcnormalized quantity plus a 
counterterm (m,- ' = rrii + Smi; g^ = g\. + 5gi). With 
these definitions, we can check that the relations ([24]) and 
(j25[) arc consistent at one-loop if we define the rcnormal- 
ized self-energy of V and <j> as 



SM V 



t v (q 2 )=j: v ( q 2 ) + 6Z v (q 2 -My) 
E (q 2 ) = E (g 2 ) + S (q 2 - to 2 ) - S to 2 . (26) 

Finally, for blobs involving a vector (V) and a scalar 
(4>) field (cf. Fig.QJ, there is a free vector index /z and 



we define % v {q) = {^) (i^v{q 2 )) = ±q»^ v {f), 
where ±q is the external momentum of </> flowing into 
(out of) the blob. As with the tensor iY 1 ^ v {q) : we equate 
the vector E^ y (q) (in this case, with no z-factor) to the 
corresponding blob with external lines (j> and V 6 . With 
these conventions, we may write the various scalar-scalar 
self-energies needed for this calculation involving both 
Higgs bosons and Goldstone bosons in diagonal or mixing 
form. For the CP-even sector: 

£<H°(<7 2 ) = ^H°(9 ) + 8 (q 2 — ?™H°) — 5 MgO 

Eh°H°(« 2 ) = E h o H o(g 2 ) + -6 Z h a H o (q 2 - m^o) 



+t;5 Z h o H a (q 2 - rri^a) - 5 ml 0R0 



2 2 \ c 2 

- m h o ) — o m h o 

2\ — *\ _/-„2\ , r r? _ / 2 2 \ c 2 



E h o( g 2 ) = Eho^ 2 ) + 6Z h o (q 2 - m£o) -Sm 2 h „. (27) 

For the CP-odd sector: 

^A°(q 2 ) = ^A°{q 2 ) + S Z A o(q 2 — m A o) — S m A o 

(q 2 ) = (q 2 ) + 8 z g° q 2 ~ S m 2 G0 
^A°G°(q 2 ) = ^A°G°(q 2 ) 



+5 Zap g° \ Q tt- ) — om 



A°G° 



(28) 



The structure of the Higgs sector beyond the leading 
order is somewhat more involved. Similarly as in the 
SM, scalar-vector mixing terms arise in the Lagrangian 
of the general 2HDM. Such contributions originate from 
the gauged kinetic terms of the Higgs doublets, 



d< = (£>„ $i) T (Z>" $i) + (Z>„ $ 2 ) T (D» $ 2 ), (29) 

where, in our conventions, the SU(2)l x U(l)y covariant 
derivative reads (using pretty standard notations), 

D lt = d lt + i^— W«T a + i^— Y B„ , (30) 

with — e being the electron charge. Upon spontaneous 
EW symmetry breaking, the following scalar-vector mix- 
ing terms are generated: 



sv 



2 s 



ir 



(v^fo + ^(jq) W+ + h.c. 



2 SyyCw 



6 These conventions are slightly different from previous compre- 
hensive calculations of Higgs boson physics presented by some 
of us in the past Q (basically they differ on the sign of the self- 
energies) and tend to be more in agreement with the conventions 
of the standard packages FormCalc and FeynHiggs 80, 81]. But 
there are some small differences; for instance, the standard func- 
tion S 'elf Energy [S — > V, Ms] of FormCalc corresponds precisely 



to our —i E 



4> V- 
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The subindex in Cq V means that this is a bare La- 



(0) 



grangian, and therefore all fields in it are bare fields 
and the parameters are bare parameters (couplings g, 
and masses rn\ ) - although we did not bother to label 
them as such in order to avoid an increasing notational 
snarl-up. 

Similarly for the VEV's V{ of the two Higgs doublets. 
Here we have to include the corresponding WF renormal- 
ization constants from (|21[) . Therefore, the decomposi- 
tion of the bare VEV's into the renormalized ones plus 
countcrterms reads as follows: 



^ (0) = Zl {vi+Svi) 
1 
2 



1 p Svi 
I 1 + - 8 + — 



(32) 



As we shall explain below in more detail, we will adopt 
the following renormalization condition for the parameter 
tan0 = «s/«i [H[8|: 



The specific structure of the various mass and WF 
counterterms becomes fixed after we choose the appro- 
priate renormalization conditions in a given subtraction 
scheme. In particular, a set of prescriptions for the 
determination of the Higgs mass countcrterms arising 
in Eq. (|28p is called for. Following the on-shell pre- 
scription, we arrange that the four renormalized Higgs 
boson masses coincide with the four input (physical) 
Higgs masses, and hence we enforce them to be the pole 
masses of the corresponding renormalized propagators 
(f2l)j). Equivalently, we declare that the real part of the 
corresponding renormalized self-energies vanishes for on- 
shell Higgs bosons. In a nutshell: 

5ReE h o (M h 2 o) = 
KeS H o (M| ) = 
5ieE A o (M| ) = 



3?eS H ± (Mh± ) = 0. 



(37) 



8 x>\ 

Vl 



8 V2 



(33) 



This condition insures that the ratio vijv\ is always 
expressed in terms of the true vacua (that is, the 
VEV's resulting from carrying out the renormalization 
of the Higgs potential) . The corresponding counterterm 
tan (3 — > tan (3+8 tan (3 associated to this definition there- 
fore reads 



8 tan j3 8v2 8v\ 1 . 

S~ = 1- o (<^*2 

tanp V2 vi 2 

= 2 (** Z ^ ~ sz <s>i) ■ 



8 Z^) 



(34) 



Proceeding in this way, we may rewrite Eq. (|31[) as Cq = 
C sv + 8C SV , with 

C sv = iM w W+ 0"G _ + h.c. - M z Z° G° (35) 

and where all the fields and parameters here are the 
renormalized ones. Similarly, the overall counterterm La- 
grangian reads 




iM w W+ d^G~ +h.c. 



iMwW+d^TT +h.c. 



M z Z° G° 
8 Z G o A o + sin 2/3^^ M z Z° d» A . 



(36) 



where we have reabsorbed the parameter v = vf + v, 
into the gauge boson masses, Mw — gv/2, Mz = 
Mw/c-w ail d their associated countcrterms. 



Incidentally, let us point out that this is a different ap- 
proach with respect to the MSSM (see e.g. [Hj]), in which 
case and owing to the supersymmetric invariance, less 
freedom is left over to settle independent renormaliza- 
tion conditions. 

However, additional conditions are needed due to the 
extended structure of the Higgs sector of the 2HDM with 
respect to the simpler SM case. In particular, the con- 
dition Q33p has been applied, but we still have to fix the 
structure of the counterterm 8 tan (3 in (|36[) , which we 
will unravel later on below. 

The scalar- vector mixing contributions at the tree-level 
in (|3"6"|) are canceled by the addition of gauge-fixing terms. 
Using the 't Hooft-Feynman gauge choice £ = 1 in the 
class of linear gauges, we have: 



C 



GF 



1 



I 2 



jr± = a M W ± ± iMw G ± 
T° = 0"Z? + M Z G° . 



(38) 



Let us, however, point out that the scalar-vector mixing 
can always take place via quantum corrections - see the 
diagrams of Fig. [1] and thereby a related set of coun- 
terterms must be introduced in order to renormalize the 
corresponding mixed self-energies. In particular, from 
Eq. dMI it follows that the A - Z° mixing counterterm 
reads 



SZ 



A°Z° 



8 Z G o A o 



1 . na 6t a n(3 
- sin 2p — 

2 tan/3 



sin 2/3 (8Z<j, 2 — SZ®^ , 



(39) 



where in the last equality we have taken 8 Z A o qq from 
Eq. (|22]l and used (|34|). Let us also derive the explicit ex- 
pression for the renormalized A — Z° mixing self-energy, 
as displayed in Fig. [TJ,. By plugging the corresponding 
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FIG. 1: Renormalized scalar- vector and scalar-scalar self- 
energies at all orders in perturbation theory. 



Feynman rules, we may cast conveniently the amplitude 
in the form 



A 



A"Z U 



M 2 

A 



q v Y, A o z o(q 2 ) 



r 



M% 



(40) 



where the structure of the renormalized A° — Z° self- 
energy encapsulates both the corresponding bare self- 
energy contribution and the related counterterm as fol- 
lows: 

Zaoz°(q 2 ) = S A oz°(9 2 ) + M Z 6Z A o z o , (41) 

with 5 Z A a z a given by Eq. Q39p . In the same vein, a sim- 
ilar expression can be derived for the mixing with the 
Goldstone boson, G°: 



A' 



M 2 



q 2 - M| 



(42) 



with £a°g° defined in (|28| . We are now ready to tackle 
the renormalization of the Higgs-boson fields. Two inde- 
pendent conditions are needed to fix d Zq, 1 , 5 Z<$ 2 . Insofar 
as we will be dealing with the neutral Higgs sector, it is 
more convenient for us to impose such conditions on the 
A boson: 



net' A0 (q 2 



3?eE AOz0 (g 2 ) 



--M 



A<> 



o, 



o, 



(43) 



(44) 



q 2 =M' 



A" 



where the short-hand notation f'(q 2 ) = df(q 2 )/dq 2 has 
been employed. The above conditions, together with 
(|37|) . ensure that the outgoing A particle is on its mass 
shell and that the residue of the renormalized propaga- 
tor (f25j) is precisely equal to 1/[1 + £ A0 (M A0 )] = 1 - see 
section |V] Besides, the second condition guarantees that 
no mixing between the A and Z° bosons will take place 
at any order in perturbation theory as long as that A 
Higgs boson is on-shcll. Actually, also the scalar-scalar 
mixing A — G° is generated through a similar class of 
one-loop interactions (recall Fig. []J>). However, one can 



see that both mixing phenomena turn out to be related 
by the Slavnov- Taylor (ST) identity 



q 2 T, A o z o(q 2 ) + M z £ A o G o(<? 2 ) 



- M lo 



0. 



(45) 



which ultimately stems from the underlying BRS sym- 
metry of the theory [H]. Notice that Eqs. (gH) and (gSJ) 
imply that the real part of the A — G° mixed self- 
energy, too, vanishes on-shell. Let us recall that a BRS 
transformation is a local Grassmann gauge transforma- 
tion in which the infinitesimal gauge parameters read 
6u> a (x) = r] a (x)5X, with rf the FP ghost fields (scalar 
fields with Fermi statistics) and where 8X is an infinitesi- 
mal Grassmann (i.e. anticommuting) parameter. In fact, 
using the notation (X) = (0| TX |0) for a Green's func- 
tion constructed out of a product of field operators X 
on which the time-ordering operator T acts, the BRS in- 
variance is expressed as 5\{X) = 0, where S\ = 8A5/5A. 
Recalling the expression for the gauge-fixing Lagrangian 
Eq. (|38[) . we consider the BRS invariance of the 2-point 
Green's function constructed with an anti-FP ghost field 
77 and a neutral CP-odd Higgs boson. We have 







(6 X (fA°)) = -(T°A°)S\ 



-(a"Z°A° + M Z G°A°)6X 



(-Z° 9" A + M z G°A°) = 0. (46) 

Here we have used the standard BRS transformation 
of the anti-FP ghost field, which is simply related to 
the gauge fixing term J 70 as follows: 5\fj° = —J 70 5X. 
Moreover, the second term of the BRS variation van- 
ishes because S\A cx A (T a /2)r] a (x)5X and we cannot have 
Green's functions with external ghost legs. 

Fourier transforming the expression (|46p . we can 
rewrite it in momentum space, 



-iq^A 



A"Z" 



Mz A 



0. 



(47) 



A" 



with A* Y denoting the renormalized (non-amputated) 
2-point Green's functions introduced in Eqs. (|4T)|) - (|4"2"j) . 
The ST-identity (|45|) is just a transcription of (|47| in 
terms of the renormalized mixed self-energies. In words, 
the identity states that an on-shell A does not mix nei- 
ther with the neutral gauge boson Z° nor with the neu- 
tral Goldstone boson G°, provided that the conditions in 
Eqs. (O and (gj) hold. A similar ST identity can be 
derived in the charged Higgs sector (see Ref. [§] ) , and al- 
though we do not need it in the context of our calculation 
we quote it for completeness. In the present notation, it 
reads 



q 2 T, H ± w ± {q 2 )±i M w S H ± G ± (q 2 ) 



q 2 = M 



0. (48) 



H± 



For an extended discussion of these identities, see 
Ref. [86| . and in particular for the off-shell regime and 
its form in generalized non- linear gauges [83| . 
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A comment on the physical definition of tan j3 is now 
in order. Recall that we have assumed that the relation 
([33]) holds between the countcrterms of the VEV's. Such 
relation says that the renormalized value of tan (3 is de- 
fined to be the ratio of the renormalized values of V2 and 
v\. This is tantamount to enforce the Higgs tadpoles to 
vanish H& 



?h {h } + $ *h° {h°} = . 



(49) 



Indeed, the tadpoles are absent at the tree-level and must 
be absent too at one-loop if v\ and vi are to be the renor- 
malized VEV's characterizing the true ground state of the 
potential at the given order of perturbation theory. 

A note of caution may be pertinent at this point: renor- 
malized quantities are not always physical quantities, al- 
though sometimes they are; for example, on-shell renor- 
malized masses are of course physical masses, and this 
is indeed what is stipulated by the fundamental on-shell 
rcnormalization conditions (|37|) . However, the renormal- 
ized VEV's Vi are not directly measurable quantities. So, 
even though the above definition of tan f3 may look quite 
natural, it is not entirely physical yet, in the sense that it 
is not tied to a physical observable at this point. In order 
to link tan /? to a genuine physical quantity we have to 
specify the process from which tan (3 could be extracted 
from the experiment. For example, in Q it is derived 
from the charged Higgs decay H — > r v T , whereas in 
other cases it is collected from other observables, e.g. 
from A a — > t + t~ [HI - see also (83[. The ultimate con- 
nection of the current definition of tan f3 with physics 
demands to take account of the contribution from the 
process-dependent quantum effects, as explained in de- 
tail e.g. in Q. Notice that the fact that this reference 
focuses on the MSSM is not relevant for this feature. All 
that said, let us stress that these process-dependent ef- 
fects are not crucially important for the present analysis 
because they do not involve the triple Higgs boson cou- 
plings under study. Therefore, since we wish to place 
ourselves in the region of parameter space where these 
couplings are most favored, the process-dependent cor- 
rections cannot noticeably alter the bulk supply of the 
significant renormalization effects on the cross-sections 
([3]) computed under these conditions. Thus, it is not 
crucial for this study of the leading triple Higgs boson 
self-interactions to deal with the details of the particu- 
lar process from which one could eventually extract the 
value of tan/3, and for this reason we will not include 
these corrections in our analysis 7 . In short, in all practi- 



7 The situation is, in the quantitative aspect, quite similar to the 
analysis of the branching ratio of the char ged Higgs decay of the 
top quark (t -> H+b) both in the 2HDM [87| and in the MSSM 
("a window to virtual SUSY" - see Ref. [1). The potential quan- 
tum effects related to the intrinsic dynamics of this decay can be 
so large that the process-dependent corrections associated to the 
definition of tan /3 pale in comparison and are, therefore, unessen- 
tial to demonstrate the very origin of the quantum corrections. 



cal respects the definition of tan f3 employed here, which 
is based on the condition ([53"]) . is well-suited to illustrate 
the leading renormalization effects produced by the Higgs 
boson self- interactions on the 2HDM cross-sections ([3]). 

Combining the renormalization conditions (|43|) and 
(|44|) . and using ([22]), ([28]) and ([39]), we are lead to 



Ke£' A o A o {M%) + sin 2 (3 SZ^ + cos 2 /3 SZ^ = (50) 

3?eS A o z o (Af 2 „) - - M z sin 2(3 (S Z* 2 - 5 Z 9l ) = . 

The WF renormalization constants for the Higgs doublets 
can now be solved explicitly as follows: 



5Z 9l = -5ReE A o(M 2 o) 
6Z$ 2 = -3?eS A o(M| ) 



1 



M z tan j3 



5fte£ A o z o(M 2 ), 



tan (3 
~M7 



5ReS A o z o(Af 2 ). (51) 



From ([3~4"]) and ([ST]) the specific form of the counterterm 
associated to our definition of tan f3 finally ensues: 



1 



Stan/3 
tan/? ~ M z sin 2/3 



3f?eS A o z o(Mlo) 



(52) 



Up to now, we have established a set of conditions 
that allows us to renormalize five out of the seven free 
parameters © of the 2HDM. These are the conditions 
([33]t . ([37]) and (j44j). Notice that the additional one ([43]) 
is a field renormalization condition and, therefore, it does 
not affect the definition of the parameters themselves - 
although it certainly affects the renormalization of the 
Green's functions and cross-sections. The two renormal- 
ized parameters that remain to be defined are the CP- 
even mixing angle a and the coupling A5. As far as a is 
concerned, we can split its bare value in the usual form 
a*- ) = a + 5a, where we define the renormalized value a 
by setting to zero the renormalized h° — H° mixing self 
energy at the scale of the Higgs mass in the final state. 
If e.g. the final state is the lightest CP-even Higgs boson 
(as in Figs. I3 1TT]) . then 



5ReSho H o(g 2 



M h 2 o) 



0, 



(53) 



whereas if it is the heavy CP-cvcn Higgs boson then M 2 
is replaced by Mg • In what follows we presume that the 
final state is the lightest CP-even state. It is not difficult 
to see that the renormalization condition ([5"3"]) fixes the 
mixing angle counterterm 8a in terms of the mixing mass 
counterterm (5to 2 i0h o in the expression of the renormalized 

mixed self-energy £ n 0H° ~~ see Eq. ([2"T]) . Obviously, ([55]) 
implies 

,5m^o H o = 3?eS h o H o (M^ ) + - 5 Z h o H o (M$, - M^ ) , (54) 

with S Zh°H° = sin 2a{5 tan (3/ tan (3). On the other hand, 
one may show from the renormalized Higgs potential that 



8a 



M 2 - M 2 



(55) 
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Actually, the condition (|53|) is bound to the fact that the 
external Higgs bosons (h°,H°) are to be on-shell to all 
orders in perturbation theory. To better assess the signif- 
icance of our renormalization conditions in the CP-even 
Higgs sector, we note that the quantum corrections in 
this sector induce the following non-trivial inverse prop- 
agator matrix: 



P 



h°H' 



S h oho 5] h o H o 

,2 n r2 , s 



p z - M£ + S h o h o 



(56) 



The physical Higgs boson masses 8 ensue from the real 



part of the roots of det A hf , H o 
equation 



0, i.e. of the algebraic 



P 



Mi 



H" 



Sh°h° ) ( P 2 - M h o + ShOh<) ) - S h o H o = . 



(57) 



The above equation makes clear that both conditions 
(|37)l and ((53)) are needed in order to settle the correct 
on-shell properties for the h° boson. Similarly, we have 
Eq. ([37]) and 5ie t h o R o(q 2 = M| ) = in the case of H°. 

The set of renormalization conditions that we have in- 
troduced so far enable us to anchor the full collection of 
countcrterms within the Higgs sector, in particular the 
mass counterterms i5mg ,fa^ 0G0 that remain unsettled 

in Eq. (|28p . In terms of the tadpole counterterms and 
6 tan j3 (all of them already fixed in this renormalization 
scheme), we may write their explicit expressions as fol- 
lows: 



8 ttu-.q 



2 swMyy 



[sin(/3 — a) St b o + cos(/3 — a) St 



H° 



S m&oao 



(58) 

[sin(/3 — a) <5t H o — cos(/3 — a) St^o] 

(59) 



2 swMw 

— sin (3 cos (3 M%p 8 tan /3 . 



The only parameter left is A5. However, this parameter is 
not involved at the tree-level in any of the cross-sections 
under study, which are entirely dependent on just the 
electroweak gauge boson couplings at the lowest order. 
Therefore, there cannot be any one-loop UV divergent 
quantity that needs to be absorbed into the renormaliza- 
tion of the parameter A5. In other words, this parameter 
is not renormalized at one-loop for the processes under 
consideration. It goes without saying that the situation 
would be different at higher orders, e.g. at 2-loop order 
and beyond. But we need not go that far to illustrate 
the main message of this paper and its potential impli- 
cations. In Section VII, we shall briefly return to this 
point. 



Let us recall that all Higgs boson masses are input parameters 
within the general 2HDM, unlike the MSSM case. 



V. NEUTRAL HIGGS-BOSON PAIR 
PRODUCTION AT 1 LOOP IN THE LINEAR 
COLLIDER: THEORETICAL SETUP 

Due to its vantage point as a high precision experi- 
mental machine, and hence owing to its complementar- 
ity with the LHC, it is of cardinal importance to under- 
stand in detail the phenomenology of the Higgs sector 
in the context of a linear collider (both within the SM 
and its renormalizable extensions). Quite an effort has 
already been devoted to this goal in the literature, spe- 
cially within the MSSM, see section [H] and references 
therein. Analyses of Higgs boson production in the LHC 
within the 2HDM have been considered e.g. in (23l. |88||. 
and the results are as encouragin g a s the corresponding 
calculation for the linear colliders [261 127I |29| except that 
in the latter case the experimental environment is much 
cleaner and the real possibilities to discriminate the na- 
ture of the Higgs bosons are correspondingly higher. 

Exclusive processes in the linear collider can be spe- 
cially significant. The simplest neutral Higgs boson pro- 
cesses of this kind in a linac involve the production of 
two Higgs bosons in the final state, Eq. At the tree- 
level, they are of order O(o 1 ew ) and are mediated by Z° 
exchange. The ncxt-to-simplest ones are the triple Higgs 
boson channels ([T]) . These have been studied in [27} , and 
in section IVIII we will compare their cross-sections with 
the 2H ones under the same set of conditions. 



Notice that due to CP-conservation - or, what is more, 
due to Bose-Einstein symmetry - processes with only 
two identical neutral Higgs bosons in the final state can- 
not proceed at the tree-level (neither in the SM nor in 
any of its extensions). Specifically, the combinations 
2H = h°h°, H°H°, h°H° of CP-even Higgs bosons can 
only be generated through CP-conserving box-diagrams 
at 0(a* w ), see Fig. [2j Such diagrams render up to 
a 2H ~ 5 x 1CT 6 pb at = 1.5 TeV [27|, which would 
translate into a rate of scarcely a few events per 500 fb _1 
of integrated luminosity, thus a too handicapped rate to 
provide feasible detection signals. It is worth pointing 
out that, since all couplings appearing in the box- type 
diagrams of Fig. [2] are fixed by the gauge symmetry, the 
obtained results do not depend on whether we consider 
the SM, the MSSM or the general (unconstrained) ver- 
sion of the 2HDM. 

Being these couplings generated from the electroweak 
kinetic terms (|29j) . they are fully determined by the gauge 
symmetry, and as a consequence they are formally the 
same both in the MSSM and in the 2HDM. It means 
that, if we aim at distinguishing between the two ba- 
sic types (supcrsymmetric and non-supersymmctric) of 
Higgs boson scenarios, the study of the quantum correc- 
tions is mandatory here. Even at one-loop order, this 
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FIG. 2: CP-conserving box diagrams corresponding to the process e + e 
describes the analogous process in the MSSM. 



h°h° in the 2HDM. A similar collection of diagrams 



implies to cope with a formidable calculation. 

Even though the detection of the 2H processes © 
would signify an unmistakable sign of physics beyond the 
SM, none of them is sensible to the trilinear Higgs boson 
couplings at the leading order. It means that the order 
of magnitude of the cross-sections (1 — 30 fb) is not sig- 
nificantly enhanced with respect to, say, the MSSM ones. 
What can be remarkably affected in the general 2HDM, 
however, is the size of the radiative corrections, i.e. the 
virtual quantum effects on these processes. For the triple 
Higgs boson channels |T]) and the gauge boson fusion 
ones ©, instead, already the tree- level cross-sections 
can be greatly modified by the triple Higgs boson self- 
interactions. Indeed, one can achieve sizable production 
rates up to roughly lpb, whereas the maximum MSSM 
payback lies around 10 -3 — 10 -2 pb (see e.g. Tables 5 
and 2 of Refs. [27| and [29| . respectively). These results 
can be traced back to the intrinsically different nature of 
the 3H self-coupling in the supersymmetric model (a pure 
gauge coupling) as compared to the non-supcrsymmetric 
ones. 

The use of the two-body processes ^ as a tool to 
probe the nature of Higgs bosons in the final states is 
nevertheless possible through the analysis of the radiative 
corrections. Closely related to them are the two-body 
final states e + e~ — > Z°h (with h = h° ,H°), which have 
long been known to be complementary to the former in 
the MSSM [28|. We will not consider the latter in the 
present analysis [89| . 

The theoretical setup for the one-loop calculation 
starts from the basic bare interaction vertices extracted 
from Eq. 



A° h° Z° 



cos(/3 — a) 
2 sw cw 



Z°A°^h° 



C 



A« H° Z° 



esin(^-q) ^ mh0 
2 sw cw 



(60) 



(61) 



In Fourier space, they read: 

e cos(/3 — a) 



A°h°Z° 



r A o H oz° 



2 sw cw 



(p-pY 



e sin(;3 — a) 
2 sw cw 



(p-P 



(62) 



(63) 



p and p 1 being the 4-momenta of the CP-odd and even 
Higgs boson pointing outwards, respectively. 

Indeed, one-loop quantum corrections to the e + c~ — > 
A°h° amplitude are driven by a large set of Feynman di- 
agrams involving the following subsets of contributions: 
i) the self-energy corrections to the Z° propagator and 
the 7 — Z° mixing propagator (cf. Figs. [3]|4j) ; ii) the 
vertex corrections to the A h° Z° interaction (cf. Figs. [5]- 
|BJ); iii) the loop-induced 7h°A° interaction (cf. Fig. [TJ ; 
iv) the vertex corrections to the e + e~Z° interaction (cf. 
Fig. H}; v) the box-type diagrams (cf. Fig. [9]); vi) the fi- 
nite wave- function rcnormalization of the external Higgs- 
boson legs (cf. Fig. [TU)) : and, finally, the countcrtcrm 
diagrams (cf. Fig. Hip. Of course an equivalent set of di- 
agrams would be needed to describe the complementary 
process e + e~ — > A°H°, but we refrain from displaying 
them explicitly. 

One might expect that the cross-sections under study 
could be sensitive to the model differences between type-I 
and type-II 2HDM's through the Yukawa coupling cor- 
rections (cf. e.g. the first two rows of diagrams in 
Fig. [5]). However, the leading effects are concentrated 
on the triple Higgs boson self-interactions, and being the 
latter identical for type-I and type-II, the model differ- 
ences are tiny in the relevant regions of the parameter 
space. 

Let us also mention that for our current purposes 
we can safely discard the pure QED corrections to the 
e + e~Z° vertex. These QED corrections and the pure 
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FIG. 3: Set of Feynman diagrams contributing to the e + e — > A°h° process at the one-loop level within the 2HDM. These 
diagrams account for the vacuum polarization effects on the Z° boson propagator. 



weak ones factorize into two subsets which are separately 
UV finite and gauge- invariant. Moreover, these photonic 
contributions arc fully insensitive, at the order under con- 
sideration, to the relevant 3H self-couplings on which we 
are focusing in this work. For the processes under con- 
sideration, where the Higgs bosons in the final state are 
electrically neutral, the one-loop QED effects are con- 
fined to the initial e + e _ vertex. In practice, the net 
outcome of the accompanying initial-state real photon 
bremsstrahlung is to lower the effective ccnter-of-mass 
energy, and hence to modify the shape of the cross sec- 
tions as a function of y/s. These QED effects arc the only 
source of infrared (IR) divergences in our calculation, and 
since we do not consider them the obtained scattering 
amplitude is IR finite. Finally, let us also emphasize that 
its inclusion would not change significantly the overall 



size of the radiative corrections, as in the regions of the 
parameter space where the 3H self-couplings arc domi- 
nant the QED effects are comparatively negligible (as we 
have checked explicitly). In short, they are unessential 
at this stage to test the presence of the new dynamical 
features triggered by the 2HDM in the processes under 
consideration. 

The UV divergent contributions arising from the set of 
Feynman diagrams that we have just described can be ab- 
sorbed by means of the associated counterterm diagrams 
of Fig. 1111 In particular, we need to include correspond- 
ing counterterms to renormalize the A h° Z° and A°H°j4 
basic interaction vertices at the one-loop level. The ana- 
lytical expressions ensue from Eqs. (|60H61[) after splitting 
parameters and fields in the usual way into renormalized 
ones plus counterterms. The final result reads as follows: 
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FIG. 4: Set of Feynman diagrams contributing to the e e — > A h process at the one-loop level within the 2HDM. They 
account for the vacuum polarization effects on the Z° — 7 mixing propagator. 
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(64) 



When writing down the above counterterms, one has 
to keep track explicitly of the UV divergences stem- 
ming from the Higgs-Higgs (h° — H°), Higgs-Goldstone 
(A°-G°) and Higgs-vector (A°-Z°) mixing effects which 
arise at the quantum level. By the same token, and de- 
spite there are no 7h°A / 7rI°A° tree-level couplings, we 
must also introduce appropriate counterterms in order to 
dispose of the associated UV-divergences that appear on 
account of the 7 — Z° mixing at the quantum level. The 
relevant terms are 



SC 



7 h°A 



cos(/3 — a) 8 Z 7 o 1 
2 siy c-w 2 



A 1 



(65) 



and 



sin(/3 



2 S\Y Cw 



a ) 5Z J^L A ^>^. (66) 



Likewise, we must introduce the corresponding countert- 
erms for the c + c~Z° vertex, as well as for the Z° — Z° 
and Z° — 7 self energies, for which we import the same 
expressions that hold in the standard on-shcll scheme for 
the SM, see e.g. [13, IzE], but of course including in their 
calculation the loop contributions from the 2HDM Higgs 
bosons (cf. Figs. [3] and [4}. Let us denote by .M^ the 
overall set of one- loop contributions to e + e~ — > h°A°, and 
by 8 Ai^ the associated counterterm amplitude. As we 
have discussed previously in terms of Feynman diagrams, 
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FIG. 5: Set of Feynman diagrams contributing to the e + e —¥ A h process at the one-loop level within the 2HDM. These 
diagrams describe part of the quantum corrections to the A h° Z° interaction vertex. 



we can split the one-loop amplitudes into the following 
(UV-finite) subsets; 



ijone—loop 
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and similarly for the case of the heavy Higgs boson, 
e + e~ — > A°H°. In the above equations, we under- 
stand that each amplitude is supplemented with the 
corresponding counterterm. The complete expression 
M one-ioo P = M (\) + § M (i) is frcc of ijy divergences, 



as we have explicitly checked. Moreover, at this point of 
the discussion we are keeping track explicitly of the X3H 
factors, in order to highlight those contributions which 
can be enhanced by the 3H self-couplings. For example, 
among the interference terms between the tree-level am- 
plitude and the one-loop amplitude (|67| we have the class 
of terms 0(a\ w ) and also the class of terms 0(a 2 ew X 2 H ), 
both classes of the same order in perturbation theory, 
although the latter arc expected to be comparatively en- 
hanced. We will check it numerically in the next section. 

Taking into account the above expressions, the overall 
(one-loop corrected) scattering amplitude for the process 
e + e~ — > A°h° reads as follows: 



■Me+e--KA°h° — \/ ^hS> -M e +e-->-AOh° 



M 



+5M {1 } 



>A°h° 



(1) 

e+e-->-A°h° 

(68) 



In the first term on the r.h.s. of this equation, we have 
appended a finite WF constant Z h o for the Higgs bo- 
son external leg. Actually, only the h° field receives a 
non-trivial contribution of this sort after we have chosen 
Za° = 1- Indeed, while we are entitled to do so for the 
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FIG. 6: Set of Feynman diagrams contributing to the e + e — > A°h° process at the one-loop level within the 2HDM. These 
diagrams describe the second part of the quantum corrections to the A h° Z° interaction vertex. 
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CP-odd field on account of the field renormalization con- 
dition (HI]), this is not possible for h°. The reason is that 
we use a WF renormalization constant for each Higgs 
doublet, not for each field, see Eq. (|2ip . and thus the re- 
lations (f22|) give no further room to normalize to one the 
residue of the propagator for the other fields. The non- 
trivial finite renormalization contribution Z h o ^ 1 for the 
other external Higgs boson in this process, h , must be 
computed explicitly, the result being 



1 



l + 5ReE' h0 (g 2 ) 



So H o(9 2 ) 



q2_ M 2 +KeSH0(9 2) 



(69) 

J 9 2 =M 2 Q 



Retaining from this expression only the 0(a ew X 2 H ) con- 
tributions to the scattering amplitude, we are left with 



>A°h° 



1 - 



»eS;„(M h 2 „) 



^SU>AOhO + 0(0: 



(70) 



Eq. (|72|) is obtained from (|28|) after substituting in it the 
explicit form of the counterterm ([59")) . 

Notice that since Eq. (|4"4"]l decrees that the real part 
of the function S A o z o(M^ ) must vanish, in the above 
amplitude (|72[) the only non- vanishing contribution shall 
stem from the imaginary part, ^mTiAPQo(M% ). A simi- 
lar expression may be derived for the H°A° channel: 



iM 



A°Z n ,A°G° 



e cos(,3 - a) ^M>g°{ m %>) 



P„. (73) 



The complete e + e — > A°h° amplitude at order 0{a\ w ) 
can finally be rearranged in the following way: 



A^c+o-^hOA 

1 - - 



M {0) x 

■ /V/1 e+e--)-h A X 



tan (/3 — a)- 



(i) 

e+e" 



► h°A° 



+ 5M {1 } 



>h°A°- 



(74) 



At the same time, our definition of the renormalizcd 
CP-even mixing angle implies that the h° — H° contribu- 
tion to the one-loop amplitude (|67p vanishes identically: 



M 



h"H" 



>A"ho 



A°h° 
2 



3?eS h o H o(M h 2 )X 



(0) 

c + c 



>A°h° 



Af h 2 -Af2 o+3?eEHO ( M 2 o ) 



(71) 



In turn, the finite correction associated to the on-shell A 
leg can be expressed in terms of the renormalized mixing 
self-energies in the guise 
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M z sin(/3 - 

Sy/CW 

e sin(/3 
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e sin(/3 - a) s a°g°(^a°) 



(72) 



where in the last step we have made use of the Slavnov- 
Taylor identity (|45[) . The 4-momentum p^ stands for 
an outgoing A boson. Finally, the expression for the 
renormalized mixing self-energy A — G° appearing in 



where the tree-level amplitude takes the explicit form: 



iM 



(0) 



2 a em 7r cos(/3 — a) 



sy/Cw ys - M z ) 
xv(pi,r)i) (# 2 - $ x ) (g L Pl + 9R Pr) u(p 2 ,m) -(75) 

In the above expression, the notation reads as follows: 
e + , with 4-momentum pi and helicity rji ; e~, with 4- 
momentum p2 and helicity 772; A , with 4- momentum 
ki; and h°, with 4-momentum k2- We have also intro- 
duced the left and right-handed weak couplings of the 
Z° boson to the electron, = ( — 1/2 + s 2 v )/cwSw- l 
9r = sy//cw and the left and right-handed projectors 
Pl.r = (1/2)(1 T 75)- Let us also mention that we have 
neglected the contribution of the Z°-width, as it is com- 
pletely irrelevant at the working energies of the linear 
colliders (s ^> M 2 ). Finally, the total cross section is 
obtained after squaring the matrix element, performing 
an averaged sum over the polarizations of the colliding 
e + e _ beams and integrating over the scattering angle: 



a(e + e 



b°A°) 



\i{s,M* ,Ml ) 
a(cost') — - x 



IM, 



32 7TS 2 
I 2 



(76) 



where we have introduced the standard definition of the 
Kahlen function X(x, y, x) = x 2 +y 2 +z 2 — 2xy—2xz~2yz. 

Before closing this section, let us stress once more that 
a similar set of formulae hold for the companion process 
e + e~ — > A°H°, which can be related to the equations 
above through the correspondence cos(/?— a) — > — sin(/3— 
a) and M^o — > Mho. 
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FIG. 7: Set of Feynman diagrams contributing to the e e — > A h process at the one-loop level within the 2HDM. These 
diagrams describe the loop-induced 7h°A° interaction vertex. 



e e e 




FIG. 8: Set of Feynman diagrams contributing to the e + e — > A°h° process at the one-loop level within the 2HDM. These 
diagrams describe the quantum corrections to the e + e~Z° interaction vertex. 
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h°\ 



FIG. 9: Set of Feynman diagrams contributing to the e + e 
diagrams describe the box-type quantum corrections. 



AV process at the one-loop level within the 2HDM. These 
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FIG. 10: Set of Feynman diagrams contributing to the e + e — > A h process at the one-loop level within the 2HDM. These 
diagrams describe the finite WF corrections to the Higgs-boson external legs. 




VI. NEUTRAL HIGGS-PAIR PRODUCTION AT 
1 LOOP: NUMERICAL ANALYSIS 

A. Computational setup 

In this section, we describe in detail the results that 
are obtained from the numerical analysis of the processes 
e+e~ -> h°A /H°A at the one-loop level. The basic 
quantities of interest are two: i) the predicted cross sec- 
tion ([76]) at the Born-level c^ ' and at one-loop ct' 0+1 ^; 
and ii) the relative size of the one-loop radiative correc- 
tions, which we track through the parameter 



In practice, the computation has been performed with 
the help of the standard algebraic and numerical packages 
FeynArts, FormCalc and Looptools (80l | for the generation 
of the Feynman diagrams, the analytical calculation and 
simplification of the scattering amplitudes as well as the 
numerical evaluation of the cross section. Aiming at a 
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TABLE III: Range of values for the A°h°Z° and A°H°Z° 
tree-level couplings - see equations (|62[) and (|63|) — that are 
used throughout the numerical analysis . 

wider survey of the regions of phenomenological interest 
in the 2HDM parameter space, we will explore the rele- 
vant range of values for both cos(/3 — a) and sin(/3 — a) 
by picking up the representative configurations quoted 
in Table Mil In particular, in the case a — (3 — tt/2 
the h°h°h° coupling takes on the SM form, as warned 
before. Furthermore, the case a = 7r/2 corresponds to 



one of the so-called fermiophobic scenarios; namely, for 
type I 2HDM, the lightest CP-even Higgs couples to all 
fermions as cos a/ sin/3 times the SM coupling (cf. Table 
Mlj) ; therefore, if a = ±ir/2 the lightest Higgs decouples 
from all the fermions 9 . Full fcrmiophobia is of course not 
possible for typc-II models (although it is partly possi- 
ble). In particular, for the MSSM this phenomenon is 
impossible and its detection would be a direct signal of 
non-SUSY physics. Let us clarify, however, that per se 
the fermiophobic scenarios are not particularly relevant 
for the current analysis because the Yukawa couplings 
do not play a central role here (as we will discuss below) . 
However, our cross-sections are indeed sensitive to the 
particular value of a that we select, and in this sense the 
two fiducial choices a = f3 and a = tt/2 represent two 
distinctive regimes for our numerical analysis, as we will 
see. 

For practical reasons, throughout our calculation we 
will be interested in projecting out the large radiative 
corrections from regions of the parameter space for which 
the predicted production rates at the tree-level are al- 
ready sizable. This means that we shall mainly dwell on 
regimes nearby the optimal value of the tree-level cou- 
pling for the IAP channel ((3 ~ a) and for the H°A one 
(f3 ~ tt/2 + a). 

In Table IIVI we quote the different Higgs-boson mass 
spectra that we shall consider in our numerical analy- 
sis. These mass sets have been designed to cover a wide 
span of phenomcnologically motivated regimes. Because 
of their mass splittings, Sets I-III can only be realized 
within a general 2HDM, whereas Sets IV- VI reflect char- 
acteristic benchmark scenarios of Higgs boson masses 
within the MSSM; they have been generated through 



9 In the two additional realizations of the 2HDM [49ll (i.e. those be- 
yond the canonical ones indicated in Tablc llTT)l - usually denoted 
as type-F and type-IP models-, the leptons couple differently 
from the quarks of the same third weak-isospin component, and 
as a result not even the type F models can exhibit full fcrmio- 
phobia. 
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the MSSM parameter inputs specified in Table [V] Note 
that Sets I to III (similarly, Sets IV to VI) are organized 
through increasing values of the Higgs boson masses. 
Furthermore, let us also highlight that, due to the pres- 
ence of a light charged Higgs boson, Sets I-II and IV-V 
are adequate only for type-I 2HDM, whereas Sets III and 
VI are valid either for type-I or type-II 2HDM. 





M h0 [GeV] 


M H o [GeV] 


M A0 [GeV] 


M H ± [GeV] 


Set I 


100 


150 


140 


120 


Set II 


130 


150 


200 


160 


Set III 


150 


200 


260 


300 


Set IV 


95 


205 


200 


215 


Set V 


115 


220 


220 


235 


Set VI 


130 


285 


285 


300 



TABLE IV: Choices of Higgs masses that used throughout our 
computation. Due to the value of M H ± , Sets I-II and IV-V 
are only suitable for type-I 2HDM, whilst Sets III and VI can 
describe either type-I or type-II 2HDM. Moreover, Sets IV- 
VI have been specially devised to reproduce the characteristic 
splitting of the Higgs masses within the MSSM at one-loop. 
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Set VI 


tan /3 


3.7 


20.0 


20.0 


M A o 


200 


220 


285 


Msusy 


300 


650 


800 




300 


300 


300 


M 2 


300 


200 


300 


X t = A t - fi/ tan/3 


-300 


-300 


-1100 



TABLE V: Choices of MSSM parameters which give rise to 
the Higgs mass sets IV- VI in Table 1 1 VI thus mimicking the 
characteristic one-loop mass splittings within the Higgs sector 
of the MSSM. The numerical value for these masses has been 
obtained with the aid of the program FemHiggs by taking the 
full set of EW corrections at one- loop [81j . Universal trilinear 
couplings [A t = Ah = At) and GUT relations for the gaugino 
soft-SUSY breaking mass terms are assumed throughout. 

In turn, the value of A5 is severely restrained by the 
theoretical constraints stemming from the tree-level uni- 
tarity and vacuum stability, i.e. from equations (|17p and 
(HUl-lgni)- The bottom line is that A 5 > is strongly 
disfavored, and the permissible values lie mostly in the 
range 

-H£A 5 <0. (78) 

The value of tan/3 that is preferred so as to optimize 
the cross-sections in the above range is tan (3 ~ 1 . We 
will carefully derive these constraints in what follows and 
shall examine their implications on the numerical analy- 
sis of the cross-sections under study. 

B. e+e" -> h°A° 

A fairly generous overview of the numerical results for 
the Higgs-pair production through c + c~ — » h°A° at one- 



loop is provided systematically in Figs [T2l - [T4l and in Ta- 
ble EH In the latter, we display the predicted values for 
the total cross-section er(c+e~ — > h°A°), together with 
the relative size of the associated quantum corrections 
([77]) for each of the Higgs-boson mass sets in Table IIVI 
The results are obtained for fixed tan/3 = 1 and differ- 
ent values of a indicated in the table. The fiducial ILC 
center-of-mass energy is taken to be either yfs = 500 GeV 
or yfs = 1 TeV. For each mass set, A5 is allowed to take 
its largest (negative) attainable value indicated in Ta- 
ble I VIII As we will confirm throughout the numerical 
analysis, this prescription insures maximally enhanced 
quantum corrections. 

The corresponding loop-corrected cross sections 
CT (o+i)^ e + e - _^ h°A°) stay in the approximate range 
2 — 30 fb for = 500 GeV - this would entail up to 
barely 10 3 — 10 4 events in the standard segment 500 fb -1 
of integrated luminosity. At yG = 1 TeV, however, the 
predicted yields deplete down to 0(1) fb in most cases, 
although this would still give a turnover of a few hundred 
events at the end of the luminosity shift. 

As for the radiative corrections themselves, they are 
also explicitly quoted in Table I VII One can see that 
they can be remarkably large, being either positive (at 
"low" energies, viz. for ^fs ~ 500 GeV) or negative (for 
y/s ~ 1 TeV and above). In the most favorable instances, 
such corrections can boost the cross-section value up to 
6 r ~ 100 % at y/s = 500 GeV, whereas the corresponding 
suppression in the high energy range (y/s = 1 TeV) can 
attain 5 r ~ —80% owing to a severe destructive interfer- 
ence between the tree-level and the one-loop amplitudes. 

In Figs [T2] - HU we explore in more detail the beha- 
vior of the radiative corrections and their interplay with 
the theoretical constraints associated to the perturbative 
unitarity and vacuum stability conditions. Although the 
range tan (3 1 is usually the preferred one from the the- 
oretical point of view, we entertain the possibility that 
tan j3 can be slightly below 1 in order to better assess the 
behavior around this value. In these plots, we depict 5 r 
(in %) for two choices of energies, yfs = 500 GeV (left 
panels) and \fs = 1 TeV (right panels). For each of these 
figures we show independent contour plots correspond- 
ing respectively to the Sets II-IV of Higgs boson masses 
(cf. Table ITV]) in the tan/3 — A5 plane. The first of these 
sets involves an assorted spectrum of Higgs boson masses, 
which is impossible to reproduce in the MSSM, whereas 
the second set closely mimics a typical MSSM-likc Higgs 
mass spectrum, in which we recognize the characteris- 
tic (approximate) degeneracy of the heavy Higgs boson 
masses: M H o ~ M^p ~ M H ± . We find that these scenar- 
ios are well representative of the phenomenological trends 
shown by all Higgs mass spectra under analysis. By set- 
ting a = (3 we optimize the tree-level A h° Z° coupling 
and in this way we insure a sizable value of the lowest or- 
der cross-section, which is of course a good starting point 
for having a chance to eventually measure quantum cor- 
rections on it. 

From these plots, it is eye-catching that the unitar- 
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FIG. 12: Contour plots of the quantum corrections S r (in %) - defined in Eq. (|77fl - to the e + e~ — > A°h° cross-section as a 
function of tan (3 and A5, for Set II of Higgs boson masses (cf. Table HV|. a = /3 and ^/s — 500 GeV (left panel), ~/s — 1000 GeV 
(right panel). The shaded area at the top stands for the region excluded by the vacuum stability bounds, whereas the disjoint 
shaded area occupying most of the central and lower region signals the domain excluded by the tree-level unitarity bounds. 




FIG. 13: As in Fig.rfJ] but for the Set III of Higgs boson masses (cf. Table [TVJ), a = /3 and yfs = 500 GeV (left panel), 
x/s = 1000 GeV (right panel). 



ity constraints are the most restrictive ones, and tend to 
disfavor large values of tan /? > 1 or values tan /? <C 1 . 
This is natural because the triple and quartic Higgs self- 
couplings rocket fast with large and small tan /?, see 
Eq. (|12p and Table HU At the same time, these unitar- 
ity constraints place a lower bound on another highly 
sensitive parameter of this study, viz. A5, whereby A5 
cannot be smaller than —11 or thereabouts (in other 
words, IA5I < 11). Besides, this parameter is sharply 
stopped from above near A5 > as a direct consequence 
of the vacuum stability condition (cf. the shaded band 



in the upper part of the plots). Thus, the combined set 
of bounds build up a characteristic physical domain in 
the parameter space of Figs [T2] - HU which pretty much 
looks like a deep valley flanked with sharp cliffs on each 
side, and centered at tan (3=1. As a result, large depar- 
tures from this central value are incompatible with large 
values of | Ag | . In particular, the range from tan/3 = 3 
onwards is circumscribed to a narrow-edged neck around 
A5 ~ 0. In addition, the heavier the Higgs mass spec- 
trum, the more severe the unitarity bounds are. Let us 
also remark that all these constraints over |Ag| exhibit 
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FIG. 14: As in Fig.[H but for the Set IV of Higgs boson masses (cf. Table HV) . a = /? and = 500 GeV (left panel), 
v/s = 1000 GeV (right panel). 





y/s = 500 GeV 


V'i = 1 TeV 




a = P 


a = /3- tt/3 


a — P — 7r/6 


a = tt/2 


a = 13 


a = P-tt/3 


01 = P — 7r/6 


a = n/2 


Set I 


Umax [ft>] 


34.13 


13.12 


24.96 


18.93 


2.89 


0.99 


1.64 


0.99 


8 r [%\ 


3.27 


58.83 


0.70 


14.57 


-72.99 


-62.69 


-79.58 


-81.47 


Set II 


Umax [fb] 


26.71 


7.34 


20.05 


13.10 


4.08 


0.85 


2.70 


1.56 


S r [%\ 


31.32 


44.43 


31.42 


28.81 


-58.42 


-65.14 


-63.28 


-68.11 


Set III 


Umax [fb] 


11.63 


3.60 


9.08 


6.36 


6.11 


1.39 


4.22 


2.86 


6V [%\ 


35.17 


67.59 


40.68 


47.86 


-30.16 


-36.36 


-35.62 


-34.58 


Set IV 


Umax [fb] 


27.44 


12.12 


18.37 


15.41 


5.99 


1.19 


2.65 


1.67 




12.86 


99.42 


0.76 


26.81 


-40.53 


-52.78 


-64.96 


-66.81 


Set V 


Umax [fb] 


23.03 


9.17 


16.90 


13.21 


6.35 


1.35 


3.48 


2.26 


5 r [%J 


21.87 


94.22 


19.25 


39.80 


-34.22 


-44.05 


-51.98 


-53.21 


Set VI 


ffmai [fb] 


8.78 


2.33 


6.86 


4.63 


7.33 


1.63 


4.78 


3.34 


S r [%\ 


19.17 


26.21 


24.13 


25.58 


-14.08 


-23.71 


-25.28 


-21.70 



TABLE VI: Maximum total cross section cr (0+1) (e + e" -> A°h°) at y/s = 500 GeV and y/s = 1 TeV, together with the relative 
size of the radiative corrections, defined by S r in Eq. (|77[) , for the different sets of Higgs bosons masses quoted in Table IWl 
The results are obtained at fixed tan/3 = 1 and different values of a, with As at its largest attainable value, cf. Table IVTll . 



a mild dependence on the actual value of the CP-even 
mixing angle a in the regimes considered here, whereas 
our cross-sections do vary significantly with it. 

In short, from the analysis of Figs. [TW+fl two basic 
regimes of phenomenological interest can be sorted out: 
1) On the one hand, we find scenarios in which A5 < 
and where this coupling can stay moderately large in ab- 
solute value (viz. | A5 1 ~ 5 — 10) while tan/3 ~ 1; this 
would correspond to the relatively narrow allowed stretch 
on the left-hand-side of these plots. In such configura- 
tions, a subset of 3H self-couplings are remarkably sin- 
gled out - the more negative is A5 , the greater is the cou- 
pling strength. Such enhancement, which is transferred 
to the e + e~ — > A°h° amplitude through a set of Higgs 
boson- mediated one-loop vertex diagrams (see Fig. [5]), 



translates into sizable quantum corrections which become 
more vigorous with growing | A5 1 . For instance, from 
Fig.[l2]we see that around tan j3 ~ 1, 6 r reaches ~ +10 % 
at A5 = —6, whilst it becomes +30% for A5 = —10. 
Furthermore, we encounter that positive radiative cor- 
rections as big as 5 r ~ +40% at \fs = 500 GeV may 
switch drastically into large negative effects of the same 
order ~ —40% (hence suppressing significantly the cross- 
section) when the center-of-mass energy is increased up 
to s/s = ITcV - sec e.g. the right panel of Fig. [12] Obvi- 
ously, positive radiative corrections are preferred in prac- 
tice because they invigorate the physical signal. There- 
fore, the startup energy y/s = 500 GeV of the ILC seems 
to be the ideal regime to probe these particular quantum 
effects rather than moving to higher energies; 
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Ecm [GeV] Ecm [GeV] 

FIG. 15: Total cross section a(e + e~ — > A°h°) (in fb) at the tree-level and at one-loop (upper panels) and relative one-loop 
correction & r (in %) - see Eq. ((77)) - (lower panels) as a function of yfs for Set II of Higgs boson masses, cf. Table [TVl Shown 
are the results obtained within three different values of A5, at tan /3 = 1 and for a — /3 (left) and a — n/2 (right) - the latter 
defining the so-called fermiophobic limit for the h° boson (for type I 2HDM) . 





A 5 


Set I 


-11 


Set II 


-10 


Set III 


-8 


Set IV 


-9 


Set V 


-9 


Set VI 


-6 



TABLE VII: Approximate maximum (negative) value of the 
parameter A5 which can be attained for each of the mass sets 
in Table HVI in accordance with the stringent constraints aris- 
ing from unitarity and vacuum stability bounds. As both sets 
of constraints depend explicitly on the input Higgs masses, we 
find a different maximum Ag j for each of the mass sets. 



2) On the other hand, a very different behavior occurs 
in the neck-shaped region of the contour plots at moder- 
ate tan/3 > 1, as it is apparent in Figs. [T2TfFH Here the 
enhancement capabilities of the leading 3H self-couplings 
are rather meager with respect to the aforementioned 
large- 1 A5 1 scenario. The reason is that for tan/3 > 3 the 
coupling I A5 1 is nailed down to take very low values. 

In the absence of a significant yield from the Higgs 
boson-mediated loop corrections, we might expect alter- 



native sources of quantum effects to pop up. The first 
ones coming to mind are of course those originating in 
the Higgs-fermion Yukawa sector. However, it turns out 
that the Yukawa interactions are not particularly effi- 
cient in the regions explored in Figs. [121X11 Indeed, from 
Table IIIII it is apparent that none of the Higgs boson 
couplings to fermions is spurred on for tan /? > 1 (in any 
of the two canonical 2HDM's). Only the regions with 
tan P < 1 - which are disfavored by the general theoret- 
ical expectations - and cos a ~ 1 could contribute here 
mainly via the Higgs-top interaction, whose strength is 
h°t t ~ m t cos a / sin j3 regardless of considering type-I or 
type-II models. Additionally, the Higgs-bottom Yukawa 
coupling (~ mfe tan/3 within type-II 2HDM), could fur- 
nish a competitive source of enhancement, as it does e.g. 
in the MSSM. In the present case, however, no trace is 
left of such effect owing to our main focus on the range 
tan/3 > 1, although this situation could change in the 
region tan/3 < 1 (see below). 

Closely related with the previous comment, we should 
also stress the fact that there are no outstanding dif- 
ferences between the obtained results from type-I and 
type-II 2HDM's in the context of this analysis. Indeed, 
to start with the Higgs-quark Yukawa interactions do not 
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FIG. 16: Total cross section <r(e + e - — > A°h°) (upper panels) and relative one-loop correction S r (lower panels) as a function 
of \fs for Set III of Higgs boson masses, cf. Table HVl Shown are the results obtained within three different values of A5, at 
tan/3 = 1 and for a — j3 (left) and a — tt/2 (right). 
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e+e"Z u [%] 
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Set I 


0.82 


-10.09 


16.56 


0.11 


7.60 


-7.00 


Set II 


0.77 


-9.84 


44.95 


0.04 


7.64 


-7.69 


Set III 


0.71 


-6.14 


35.15 


-0.23 


7.40 


-0.57 


Set IV 


0.70 


-7.37 


20.21 


-0.30 


7.59 


-5.79 


Set V 


0.68 


-7.12 


26.32 


-0.30 


7.62 


-3.49 


Set VI 


0.66 


-6.26 


15.06 


-0.20 


7.69 


2.64 



TABLE VIII: Relative contribution to the total cross section 8 r> i = (cr 0+l ' — a^°')/a^ a ' from each topology of one-loop diagrams 
normalized to the tree- level rate. The results are derived at fixed yfs = 500 GeV for tan/3 = 1, a — /3 and choosing the value 
of A5 that optimizes the quantum corrections for each mass regime, see Table IVIII 



enter at the tree-level in the processes ((3]) under consi- 
deration; second, the trilinear couplings arc common to 
both types of models; and third, we have already empha- 
sized that the one-loop differences that could manifest 
through the distinct form of these Yukawa couplings in 
type-I and type-II models are virtually obliterated in the 
physical region permitted by the constraints. All in all, 
perhaps the most distinctive feature between both mod- 
els boils down, in this context, to just the requirement 
that M H ± > 300 GcV for typc-II 2HDM's as a result 
of the low-energy B-mcson physics constraints, which 
tend to substantially raise the average mass of the Higgs 
spectrum for this class of models - cf. Sets III and VI 



of in Tabic IIVI The basic differences in the physics of 
e + e _ — > A°h° in both models thereby narrow down to 
mass differences between their characteristic Higgs-boson 
mass spectra 10 . 

A very similar pattern is met in the other explored 



Recently, a combined analysis of different B-mcson physics con- 
straints over the 2HDM parameter space suggests that values of 
tan /3 ~ 1 could be disfavored for charged Higgs masses M H ± too 
near to 300 GeV (and certainly below) [45| . However, the level 
of significance is not high and, moreover, our leading quantum 
corrections are basically insensitive to the charged Higgs mass. 
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FIG. 17: Total cross section <r(e + e - — > A°h°) (upper panels) and relative one-loop correction 5 r (lower panels) as a function 
of \fs for Set IV of Higgs boson masses, cf. Table [TV] Shown are the results obtained within three different values of A5, at 
tan/? = 1 and for a — j3 (left) and a — tt/2 (right). 



scenarios, the only differences being ascribable to the dis- 
tinct Higgs boson mass spectrum that each of these sce- 
narios accommodate. In particular, heavier Higgs boson 
masses may have a twofold impact: i) dynamical, in that 
they soften the influence of Higgs-mcdiatcd one-loop di- 
agrams - and hence the strength of the 3H self-couplings 
- inasmuch as these contributions become suppressed by 
inverse powers of the Higgs masses; and ii) kinematical, 
in that they shift the peak of the (y{\/s) curve, since the 
latter is correlated with the production threshold of the 
Higgs pair. 

All the above mentioned phenomenological features are 
transparent in Figs. [T51IT71 in which we explore the evo- 
lution of the cross section as a function of the center- 
of-mass energy. We include, in each plot, the tree-level 
contribution and the loop-corrected one cr^ 0+1 - ) for 
different values of A5. Worth noticing is that by testing 
the influence of A5 on the cross-sections, we are testing 
the bulk capability of the 3H self-interactions, and hence 



Therefore, a shift of M H ± slightly upwards restores the possi- 
bility of tan f3 ~ 1 at, say, 2cr without altering significantly our 
results (as we have checked.) 



their potential fingerprint on the cross-sections at the 
quantum level. In the lower panels, we track the related 
behavior of the quantum correction 5 r with ^fs for the 
same set of A5 values. The plots are generated at fixed 
tan j3 = 1 and for two different values of the tree-level 
A h° Z° coupling: namely, a = /3 (where the tree-level 
coupling is maximum) and a = 7r/2 (corresponding to 
the aforementioned fermiophobic limit, in which h° fully 
decouples from the fermionic sector for type I models). 
In Figs. 1151161 and [TTl we display the numerical results 
obtained for Sets II, III and IV of Higgs boson masses, 
respectively. A similar numerical output is obtained for 
the remaining mass sets in Table IIV1 showing no relevant 
departure from the phenomenological trend that we have 
recorded as yet. 

Let us focus e.g. on Fig. [15] for a while. The tree- 
level curve exhibits the expected behavior with yfs, 
as it scales with the s-channel Z°-boson propagator ~ 
l/(s— M|). The maximum cross section is achieved at 
~ 500 GeV. The same pattern is followed by the full loop- 
corrected cross-section, although the dependence with y/s 
is sharper. In fact, the range where S r > is much 
briefer than the one where 5 r < 0, although the former 
is characterized by very significant quantum effects; so 
much so, that they may skyrocket up to < 100 % (at the 
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FIG. 18: Total cross section <r(e + e - — > A°h°) (upper panels) and relative one-loop correction 8 r (lower panels) at fixed 
V's = 500 GeV, given as a function of M h o (left) and M A0 (right) for Set II of Higgs boson masses, cf. Table [TV] tan/3 = 1, 
a = P and three different values of A5. 



largest allowed values of | A5 1 ) around the critical energy 
domain where yfs = 500 GeV. Such effects could hardly 
be missed in the first runs of the ILC. In contrast, when 
we move from y/s = 500 GeV to y/s = 1000 GeV the re- 
sulting cross section <j(° +1 ) drops by at least 50%. The 
impact of the enhanced 3H self-couplings can be easily 
read off from the dramatic differences in the one-loop pro- 
duction rates when the value of A5 is varied in the theoret- 
ically allowed range ([78]) , Thus, when |As| is pulled down 
at fixed tan/?, the quantum effects are rapidly tamed, in 
correspondence with the fact that the 3H self-couplings 
are no longer stirred up. 

The drastic decline of the triple Higgs boson self- 
interactions in this domain cannot be counterbalanced 
by the contribution of the Yukawa couplings, as they 
are not very significant in the region of the parameter 
space where tan/3 is of order one. At this point, it comes 
to mind our former observation about the potential role 
that the region tan /3 < 1 could play. Even though it is 
not favored on theoretical grounds, we can not exclude 
the possibility that 0(0.1) < tan/3 < 1, and in fact the 
unitarity and vacuum stability bounds cannot dismiss in 
block this range when A5 is sufficiently small. Thus, for 
example, in the conditions of Fig. 1121 and under the as- 
sumption that A5 = 0, we have checked that tan/3 ~ 0.2 
is allowed by the aforesaid constraints, and in this region 



we encounter quantum effects as large as +30%. They 
are the result of the combined effort from the diagrams 
involving the top quark Yukawa coupling and those in- 
volving the trilinear Higgs boson couplings - which, in 
this case, become enhanced by a pure tan/3 < 1 effect, 
rather than by an oversized A5. 

It is also interesting to track carefully the fractional 
payoff S r = (er(° + ') — 0-(°))/cr( o ) from the various one- 
loop topologies (i — 1,2,3,...) of Fcynman diagrams con- 
tributing to the process c + c~ — ► h°A°. The correspond- 
ing results are displayed in Table IVIIII under the given 
conditions. Included there are the gauge invariant sub- 
sets of diagrams of various types, namely Z° — Z°/7 vac- 
uum polarization effects, box diagrams, the renormalized 
A h° Z° and e + e _ Z° vertex corrections, the loop-induced 
form factor 7h°A°, and finally the finite WF corrections 
to the external h° leg. We can see that the gauge bo- 
son self-energies provide tiny contributions of order 1% 
at most. The box diagrams, in their stead, are non- 
negligible (of order 10%), albeit negative. At the same 
time we detect mutually destructive effects between e.g. 
the finite WF correction and the e + e~Z° vertex. 

Most conspicuously, and above all, we revalidate 
the dominant role of the Higgs boson mediated loop- 
diagrams (triggered by 3H self-interactions) as the truly 
leading source of quantum corrections to the cross- 
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FIG. 19: Total cross section <r(e + e - — > A°h°) (upper panels) and relative one-loop correction 8 r (lower panels) at fixed 
y^i = 500 GeV, given as a function of Af h o (left) and M A0 (right) for Set III of Higgs boson masses, cf. Table ITVl tan ft = 1, 
a = f3 and three different values of A5. The dashed regions are excluded by the constraints over the Sp parameter. 



section. These dominant effects are concentrated on the 
fundamental A h° Z° vertex. This vertex, which is purely 
gauge-like at the lowest order of perturbation theory, be- 
comes drastically promoted at the one-loop level as a 
result of the dominant quantum effects sourced by the 
3H self-interactions. We may cast the overall renormal- 
ization of such vertex as follows 11 : 



r A« h» Z° Z l 



1/2 ( r 

1 A o h o z o 



,1/2 



1 A° h° Z° 



r A0 h° Z° 

r 1.3H 
1 AO h° Z° 



(79) 



where T 1 ' 311 denotes the leading subset of (Higgs- 
mediated) one-loop corrections. The latter can be en- 
capsulated in a UV-finite subset of one-loop Feynman 
diagrams and subsequently reabsorbed into a modified 
or effective A h° Z° coupling. Simple power counting and 
educated guess provides the following estimate: 



-<eff 

A° h° Z° 



^0 

A h° Z° 



A 



3H 



16ir 2 s 



f(M^o/ S ,M 2 AO / S ) (80) 



A similar discussion is of course also applicable to the comple- 
mentary A H° Z° coupling and the corresponding production 
process e + e~ 



H°A°, see the next subsection. 



where f(M? /s, M^ a /s) is a dimensionless form factor 
that accounts for the complete one-loop amplitude. For 
all its cumbersome structure, we do know at least that / 
dies out with the center-of-mass energy of the process and 
encodes the sign flip of the radiative corrections when we 
increase the center-of-mass energy from yfs = 0.5 TeV to 
y/s = 1 TeV. The prefactor 1/167T 2 stands for the trade- 
mark numerical suppression attached to one-loop form 
factors, whilst A 3 # denotes a generic 3H self-coupling. 
Noteworthy is the fact that this coupling can reach val- 
ues as high as \ 3 h/M w ~ | A 5 1/ e ~ 30 and still under 
the stringent constraints imposed by unitarity. 

It is enlightening to try to estimate the expected size of 
the maximum corrections. From Eq. (|80p we can roughly 
gauge the maximum relative contribution, e.g. at the 
startup energy y/s = 500 GeV: 



5 r 



r(0) 



<|7W(°)| 2 ) 



IA 



3H\ 



16 7T 2 S 



f(M* /a,Ml /s) 



\/s=500 



50%, (81) 



where < ... > stands for the various operations of aver- 
aging and integration of the squared amplitudes. In the 
above equation, we have neglected the square of the one- 
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FIG. 20: As in Fig. [19] but for Set IV of Higgs boson masses, cf. Table HVl 



loop amplitude (see, however, just below) and set / ~ 1, 
which is a reasonable assumption since the current value 
of the center-of-mass energy, = 500 GeV, is far from 
the range in which one-loop corrections are proved to be 
negligible. The latter result of Eq. ([8T]) falls in the right 
ballpark of the optimal values for S r that we have iden- 
tified from our exact numerical analysis, for Higgs boson 
masses of the order of a few hundred GeV and maximum 
allowed | A5 1 values (so as to push X^h close to the unitar- 
ity border). By means of the same formula we can also 
estimate the contribution to 5 r arising from the squared 
of the one-loop amplitude, which we have neglected be- 
fore. We find 



(l-MT) 

(l-MT) 



X 2 



I ^=500 



5%, (82) 



thus roughly 10% of the one-loop leading effect. It turns 
out that this is again in good agreement with the nu- 
merical calculation 12 . All in all we conclude that, as 
long as 3H couplings are sizable, quantum effects will be 
manifest - and potentially very large. 

Some more comments are in order concerning impor- 
tant features contained in Figs. [T5HT71 As we move from 



12 We will further dwell on the potential relevance of such 0(a^ w ) 
effects in Section I VIII 



Fig. [15] to [TCI hence as we transit from a relatively light 
to a heavier Higgs boson spectrum (cf . Table IIV[) , we 
encounter that the resulting cross sections, both at the 
leading and at one-loop level, are pulled down by roughly 
a factor of 2. Unsurprisingly, this is signaling the ex- 
haustion of the available phase space. By the same to- 
ken, the peak-shaped maximum of the <j(^fs) curve gets 
shifted forward by about 100 GeV, in correspondence to 
the production threshold of the Higgs pair. Moreover, 
the heavier the Higgs bosons, the more suppressed the 
Higgs-mediated loops become, and hence the lower is the 
influence of 3H self-couplings. This is precisely what we 
confirm numerically in Fig. [TS] (lower panels) , in which 
we can appreciate how 5 r reaches ~ 40% at most - in 
contrast to S r ~ 100% that was spotted for Set II (cf. 
Fig. IT5"j) . In turn, Set VI (in Fig. [T7|) embodies again a 
rather light Higgs boson spectrum which resembles that 
of Set II, even though the present choice of masses is ob- 
tained upon artificially engineering a SUSY setup. Here, 
owing to the fact that Af H ± is heavier than in Set II, 
unitarity conditions settle a more stringent upper bound 
on I A5 1 , and hence on the maximum enhancement of the 
triple neutral Higgs couplings. As a result, radiative cor- 
rections in this case are still sizable, but retreat to the 
- 40% level. 

Of great relevance is the study of the behavior of 
cr(e + e~ — > Ii°A°) as a function of the Higgs boson masses. 
We present these results conveniently in Figs. 
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y/s = 500 GeV 


,/i = 1 TeV 




a = P - tt/2 


a = p- tt/3 


a = p — 7r/6 


a = 


a = P — 7r/2 


a = P - tt/3 


a = p — 7r/6 


a = 


Set I 


fT moT IT) 


40.54 


20.33 


2.29 


8.82 


2.51 


1.30 


0.16 


0.59 


Sr [%\ 


51.83 


1.50 


-65.72 


-33.96 


-75.59 


-83.10 


-93.62 


-88.50 


Set II 


<Jmr,i- ID 


26.32 


17.53 


4.39 


10.29 


3.59 


2.52 


0.69 


1.55 


S r [%\ 


48.45 


31.82 


-1.03 


16.10 


-62.62 


-65.09 


-71.39 


-67.75 


Set III 


O 'max [fb] 


5.00 


3.28 


0.87 


1.95 


5.17 


3.98 


1.41 


2.75 


S r \%\ 


71.81 


50.22 


19.75 


34.15 


-36.80 


-35.14 


-30.77 


-32.80 


Set IV 


Cmax [fb] 


20.20 


10.78 


1.98 


3.77 


3.23 


4.32 


1.14 


2.16 


S r [%] 


102.73 


44.21 


-20.53 


-24.43 


-64.02 


-35.78 


-49.13 


-51.97 


Set V 


(Tmax [fb] 


11.42 


3.75 


0.98 


2.13 


4.03 


1.69 


1.04 


1.82 


S r [%\ 


115.61 


-5.53 


-26.13 


-19.69 


-52.56 


-73.46 


-51.19 


-57.22 


Set VI 


Umax [fb] 




below 
threshold 


4.76 


4.12 


1.44 


2.61 


Sr [%\ 




-26.91 


-15.71 


-11.69 


-19.95 



TABLE IX: Maximum total cross section cr (0+1) (e + e" -> A°H°) at y/s = 500 GeV and y/s = 1 TeV, together with the relative 
size S r of the radiative corrections, for the different sets of Higgs bosons masses quoted in Table IWl The results are obtained 
at fixed tan/3 = 1 and different values of a, with | As | at its largest attainable value, cf. Table IVlTI 



where again we superimpose both the tree-level and the 
loop-corrected cross-sections for the process e + e~ — ► 
h°A°. The dependence with M h o is indicated on the left 
panels and that with Mpp on the right panels. The nu- 
merical results are shown for the Sets II, III and IV and 
have been obtained, once more, by setting tan /? = 1 , 
a = j3, and making allowance for to take on the 
highest permissible value in each set (see Table IVII[) . 
The center-of-mass energy is settled at the fiducial value 
Vs = 500 GeV. 

The plots under consideration evince that the raise of 
the Higgs boson masses exerts a twofold effect on a, one 
kincmatical and the other dynamical. The former sim- 
ply means that owing to the reduction of the available 
phase space the cross-section obviously falls down. We 
can check it in the figures, where we see that both the 
tree-level and the loop-corrected cr^ 0+1 - ) indeed de- 
crease monotonically with growing M^o and Mao. As 
far as the quantum corrections themselves are concerned, 
the situation is a bit more subtle. Although heavier 
Higgs bosons imply larger propagator suppressions of the 
Higgs-mediated one-loop diagrams, as we have just de- 
scribed in the previous paragraph, the 3H self-couplings 
can partially offset this situation since they get invigo- 
rated in such circumstances. This counterbalance fea- 
ture is manifest in the explicit expressions for the 3H 
couplings in Table |H] - the uttermost limitation to their 
growing being the unitarity constraints. In particular, 3H 
couplings involving neutral fields (which can be shown to 
carry the dominant effect on c + e~ — > h°A°) turn out to 
grow either with M h o or Mpt>. Although both effects 
partially balance each other, we are ultimately left with 
the expected decoupling behavior in the limit of large 
Higgs boson masses, with the proviso that such decou- 
pling takes place in a softer way than that of the Born- 
level cross section. In fact, in Figs [TSH2"01 we can verify 
that the one-loop corrected cross-section remains most of 



the time on top of the tree-level result when we increase 
the Higgs boson masses, and moreover the slope of the 
former is clearly milder, particularly for M h o < 200 GeV. 
This translates into a growing trend of S r up to values of 
M^o ~ 275 GeV. Let us also mention in passing that the 
sudden spike in the <5 r (MAt>) curve at Map — 340 GeV 
(which is also barely visible for a) corresponds to the tt 
pseudo-threshold in the one-loop vertex amplitude (cf. 
the first row of diagrams in Fig. [5]). 

Analogous trends are encountered for sets III and IV 
of Higgs boson masses, the only reportable difference be- 
ing that the decoupling behavior of cr( 0+1 ) with M h o, Map 
is more pronounced, signaling that the increase of these 
masses does not boost the 3H self-couplings as efficiently 
as in the case of Set II. The shaded regions in Fig. [T5]|2TJ1 
correspond to mass ranges already excluded by experi- 
mental measurements of the Sp parameter (cf. section 

ED. 



C. e+e" -> A°H° 

On general grounds, the basic phenomenological fea- 
tures are common to the h°A° and H°A° channels, with 
the only proviso that, in the latter case, the result- 
ing production rates are slightly smaller simply because 
Mho > M h a. The final cross-sections are nevertheless 
comparable and reach up to a few dozen fb. In Ta- 
ble |IX] we collect, in a nutshell, the essential results from 
a dedicated numerical analysis of cr(e + e _ — > H°A°) for all 
sets of Higgs boson masses (cf. Table IIV[) and different 
values of the tree-level A H° Z° coupling, again at fixed 
tan/? = 1 and using the largest permissible value of |As| 
in each set (see Table IVII|) . The relative size of the quan- 
tum corrections 8 r is also included in Table HXl As in the 
previous channel, we have used fiducial center-of-mass 
energies y/s = 500, 1000 GeV respectively. The maxi- 
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FIG. 21: Total cross section a(e + e~ — > A°H°) (upper panels) and relative one-loop correction 8 r (lower panels) as a function 
of \fs for Set II of Higgs boson masses, cf. Table ITVl Shown are the results obtained within three different values of As, at 
tan/3 = 1 and for a = j3 — iv/2 (left) and a = (right) - the latter defining the so-called fermiophobic limit for the H° boson 
(for type-I 2HDM). 



mum production rates at yG = 500 GeV arc achieved 
for relatively light Higgs boson masses (Set I) and for 
a = f3 — 7r/2, in which case the tree-level A°H°Z in- 
teraction is not suppressed, see Eq. (|oTj). Within this 
setup, one is left with cr( 0+1 ) ~ 40 fb. Indeed the pro- 
duction rates lie at the O(10 fb) level in a broad range 
around a = /3 — n/2 for all the mass sets under analy- 
sis. If we nevertheless consider a = /3 — 7r/6, when the 
tree-level coupling gets suppressed by sin(/3 — a) = 1/2, 
the resulting cross sections are still at the level of a few 
fb - each femtobarn rending 500 events per 500 fb -1 of 
integrated luminosity. 

Radiative corrections thus may leave a solid imprint 
into the final predicted er( 0+1 ) of the current process. In 
particular, for a = (3 — ir/2 quantum effects rubber-stamp 
a characteristic 50—100% boost upon the tree- level cross- 
section, which could hardly be missed. Other choices of 
a amount to negative one-loop contributions in which 
A H° Z° is dramatically suppressed. For instance, in the 
case of Set I for a = (3 — tt/6, it renders 6 r ~ —65%. 
For center-of-mass energies from around i/s = ITeV on- 
wards, a powerful destructive interference is seen to be 
under way between the Born and the one-loop ampli- 
tudes. This translates into large (negative) quantum cor- 



rections which may well reach 5 r ~ —50 — 100%, thus lit- 
erally stamping out the signal and pulling it down at the 
~ fb level or even below. To better illustrate this feature, 
in Figs. l2"TH2"3"l we plot the behavior of a(c+e~ -> A°H°) 
(top panels) and S r (bottom panels) as a function of 
for the representative Sets II-IV of Higgs boson masses 
of Table IIVI and for different values of the tree- level cou- 
pling; viz. a — ft — 7r/2 (left panels), which corresponds 
to the maximum of the tree-level A H° Z° coupling, and 
sin a = (right panels), corresponding to the fermiopho- 
bic limit for H° boson within type-I 2HDM (cf. Table 
IIII[) . Again we fix tan/3 = 1 and we vary A5 within 
the range allowed by unitarity and vacuum stability re- 
strictions. As it was already observed in the analysis 
of the h°A° channel, heavier Higgs-boson spectra fur- 
nish smaller production rates. We refrain from extend- 
ing the discussion, which is entirely similar to that of 
c+c- -> h°A°. 

To our knowledge, no pre-existing study of the one- 
loop effects on the neutral Higgs boson channels §5$ can 
be found in the literature. Let us, however, briefly com- 
ment on the prospects for the charged Higgs-pair chan- 
nel, e + e~ — > H + H~. This was considered some time 
ago in references [4U - I43I ] . These old studies were devoted 
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to the computation of er(c + e — > H + H ) at one-loop in 
both the MSSM and the 2HDM. In the latter case, they 
also uncover very large quantum effects (of order 100% 
or even higher 13 ) which are attributed to enhanced 3H 
self-interactions involving charged Higgs bosons. Such 
quantum effects may lie well above the typical MSSM 
counterparts, which arc found to render up to \5 r \ ~ 20%, 
and mostly negative. Let us, however, remark that the 
c + c~ — > H + H~ channel could operate at sizable cross- 
sections levels only for type-I models since the type-II 
models require a rather heavy mass spectrum. 



We should nevertheless put a note of caution here in that the uni- 
tarity constraints that were employed in Ref. I41I443H are substan- 
tially less restrictive than those considered throughout our work. 
In this regard, an upgraded calculation of c(e+e — — > H+H — ) 
is mandatory before we can judiciously compare with the results 
that we have presented here for the neutral Higgs boson channels. 
However, due to the increase of the charged Higgs boson mass 
bounds in the last few years, the maximum cross-section rates 
should correspondingly fall down, thereby making this process 
less competitive. 



VII. DISCUSSION AND CONCLUSIONS 

In this work, we have undertaken a thorough inves- 
tigation of the pairwise production of neutral Higgs- 
bosons within the context of the general two-Higgs- 
doublet model (2HDM) at the ILC/CLIC linear colliders 
(linac) through the basic processes c + e~ — > h°A /H°A . 
To the best of our knowledge, this study is the first one 
existing in the literature which addresses these impor- 
tant processes at the one-loop level. To sum up the main 
objectives of our work, we can assert that we have con- 
centrated on the following three tasks: i) we have, on 
the one hand, singled out the most favorable regimes 
for the Higgs-pair production processes once the radia- 
tive corrections are taken into account, specifically within 
the framework of the on-shell renormalization scheme 
with an appropriate definition of tan/3 and of the CP- 
even mixing angle a at the quantum level; ii) next, 
we have carefully quantified the importance of the lead- 
ing quantum corrections, and not only so, we have pin- 
pointed also the origin and relative weigh of the vari- 
ous sources; and, finally, iii) we have correlated the most 
powerful source of quantum effects with the strength of 
the Higgs boson self-interactions (essentially the triple 
Higgs boson self-couplings X^h), to which the processes 
e + c~ — >• h A /H°A are highly sensitive through Higgs- 
mediated one-loop diagrams. With these aims in mind, 
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we have systematically searched over the 2HDM param- 
eter space by considering different Higgs boson mass sce- 
narios and different values of the tree-level A h° Z° and 
A H° Z° couplings. At the same time, we have carefully 
explored the very frontiers of the generic 2HDM param- 
eter space, specially the limits that cannot be trespassed 
by the strengths of the various 3H self-couplings. 

Among the various subtle triggers that control the size 
of these important couplings, one of them stands out with 
overwhelming supremacy over the others, and this is the 
A5 coupling. This is the only one parameter in the list 
of seven free parameters - cf. Eq. (j9)) - conforming 
the structure of the general (CP-invariant) Higgs poten- 
tial of the 2HDM that cannot be ultimately traded for 
a physical mass or a mixing angle. This parameter is, 
in addition, the only Higgs self-coupling that need not 
be renormalized at one loop in the current framework. 
Despite its special status of being apparently unfettered 
and without obvious physical connotations, one quickly 
discovers that it becomes severely restrained by subtle 
bounds connected with the quantum field theoretical con- 
sistency of the model, namely the stringent constraints 
imposed by tree-level unitarity and vacuum stability con- 
ditions. These restrictions, together with the limitations 
dictated by custodial symmetry, and of course also those 
emerging in the light of the direct collider bounds on the 
masses of generic Higgs boson models, not to mention the 



indirect phenomenological restraints (in particular, the 
upper bound on the charged Higgs mass in type-II mod- 
els ensuing from radiative B-meson decays), do indeed 
eventually determine a fairly compact region of param- 
eter space. In particular, A5 is compelled to be mostly 
negative and in the range | A5 1 < 0(10). At the same 
time, tan/3 cannot be large when |As| is also large with- 
out violating the constraints. The bottom line is that the 
value that maximizes the processes under consideration 
is a rather equitable one: tan/3 ~ 1. 

Notwithstanding the sharp limits placed upon the 
physical region of the parameter space, our explicit cal- 
culation of radiative corrections has shown - and we feel 
it is truly remarkable - that the general 2HDM models 
are still able to unleash a good deal of their "repressed" 
power. If only we could catch a smoking gun of this stu- 
pendous potential through the accurate measurements 
that a linear collider should be able to render on funda- 
mental processes like e + c~ — > h°A°/H A , we might be 
on the verge of experiencing an intense episode of Higgs 
boson physics beyond the SM. In the following, we de- 
scribe the basic results and strategies that we could follow 
to this effect. 

We have identified three main scenarios of potential 
phenomenological interest in the corresponding parame- 
ter space. In all of them we observe no remarkable de- 
pendence either on the details of the Higgs mass spec- 
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trum nor on the pattern of Higgs-fermion Yukawa cou- 
plings (either for type-I or type-II 2HDM), not even on 
the particular channel h A°/H°A under study. In what 
follows we describe each one of these scenarios and spell 
out clearly their differences, in particular we provide the 
characteristic potential size of the relative quantum cor- 
rection 6 r - cf. Eq. (|77| - typically associated to them 

(cf. Figs-nanj: 

• Scenario 1. To start with, regions of tan/3 > 1 
and large | A5 1 (specifically, A5 < 0, | A5 1 ~ 5 — 10) 
support the bulk of the radiative corrections. They 
may well reach S r = ±(50 — 100)% at the largest 
attainable values of ] A5 [ . Furthermore, these effects 
critically depend on the actual size of |As| and on 
the center-of-mass energy y/s (cf. Figs. lTMT7j) . to 
wit: 5 r is positive at "low" center-of-mass energies 
around y/s = 500±100 GeV, while it is negative for 
y/s > 600 GeV and certainly also in the uppermost 
foreseeable energy segments y/s =1 — 3 TeV, usu- 
ally reserved for an upgraded ILC and specially for 
CLIC. Thus, interestingly enough, quantum cor- 
rections turn out to drastically invigorate the cross 
section of the basic processes e + c~ — » h A°/H°A 
near the standard booting energy of the ILC, while 
they tend to suppress the two-body Higgs boson 
signal for a center-of-mass energy some 20 — 30% 
beyond this initial regime - and of course for the 
TeV range and above - as a result of construc- 
tive/destructive interference effects, respectively, 
between the tree-level and the one-loop amplitudes. 
Incidentally, we clarify that this optimal scenario 
cannot be iterated for large and positive values of 
A5 because the combined set of constraints (mainly 
the vacuum stability bounds, in this case) preclude 
most of the region A5 > 0; 

• Scenario 2. On the other hand, we find a tail of 
subleading effects within a band of the (tan/3, A5) 
subspace at moderate tan/3 = 1.5 — 5 and |As| < 2 
wherein the maximum radiative corrections turn 
out to stagnate and remain barely at the level of 
~ 10% or less. Should the parameters of the model 
inhabit this region, the practical possibilities to de- 
tect these effects would obviously become thinner 
than in the previous case. Still, 5 — 10% effects 
are not that alien to the high precision standards 
planned for a linac collider, and hence it should not 
deter us from the searching task. At this point, the 
detection strategies associated to the Higgs boson 
decay patterns could be essential (see below); 

• Scenario 3. Finally, there is the region of pa- 
rameter space where the coupling A5 is very small 
or it just vanishes. Here the quantum corrections 
can still have a chance provided tan/3 < 1, and 
indeed we have found domains of this kind preserv- 
ing all the basic constraints. While this region is 
not the most favored one, neither from the the- 
oretical point of view nor from the experimental 



data, specially if tan/3 <C 1, we should keep it in 
mind as a possibility. Even for moderately small 
tan/3, say for tan/? > 0.2, it can be the source of 
sizeable quantum effects of order 10 — 30% (and 
positive), which are larger than those in Scenario 
2 and, in some cases, even comparable to those in 
Scenario 1. It is worth noticing that, in this region, 
there is an interesting collaborative interplay be- 
tween the Higgs boson trilinear self-couplings and 
the top quark Yukawa coupling (which is the same 
in type I or type II 2HDM's). Only in this scenario 
a Yukawa coupling could play a significant quanti- 
tative role on equal footing with the trilinear Higgs 
boson self-couplings. 

Obviously, the success of the whole Higgs boson search 
endeavor will also depend on the pattern of distinctive 
signatures available to these particles in the final state. 
For instance, for M h < 2M V < 180 GeV, we should 
basically expect back-to-back pairs of highly energetic 
(roughly E ~ y/s/ 4 > 100 GeV) and collimatcd b-quark 
and/or T-lepton jets from A°,h° — > bb/r + T~ which, in 
contrast to the LHC, should not be overshadowed by 
the large QCD background which is inherent in hadronic 
machines. And of course similarly with the A°H° fi- 
nal state. This much clearer environment notwithstand- 
ing, we should be aware of the fact that other char- 
acteristic processes of linac physics, such as the four- 
jet cross-section, might be large enough to partially ob- 
scure the Higgs-pair production signal. However, the pre- 
cise analysis of the corresponding signal distributions, 
which would be required in order to completely assess 
the real experimental possibilities, is beyond the scope 
of the present study. And of course the same reason- 
ing applies to the A°,H° final state. For M h > 2M V , 
on the other hand, signatures from h°, H° — > VV and 
A — > 66/t + t~ could also play a role (depending on 
the particular choice of a and /3, cf. Table HIT) . They 
could ultimately lead to signatures with two (up to four) 
charged leptons against a bb or t + t~ pair in the final 
state (from W ± — > £ + missing energy and, specially, 
from Z — > £ + £~). Even if the branching ratios of gauge 
bosons into leptons are very small (B(V —> 11) ~ 0.03), the 
predicted number of Higgs-pair events is large enough 
so that, in the most favorable scenarios, one could col- 
lect O(10 2 ) events for these complementary signatures, 
which could be instrumental for a final tagging and iden- 
tification of the Higgs bosons. Notice that we use the 
fact that the CP-odd A "always" decays into bb or, to a 
lesser extent, to t + t - (even for tan/3 = 1) as it cannot 
couple to gauge bosons. So the patterns of signatures are 
rather characteristic and could not be missed in normal 
circumstances. But there is more than meets the eye in 
the whole strategy game, as we shall see later on. 

In the most favorable regimes (Scenario 1), the exper- 
imental chances could be spectacular, especially in the 
case of relatively light Higgs boson mass spectrum and 
a large (negative) value of As, for which the predicted 
one-loop cross sections a^ 0+1 \e + e~ h A°/H A ) at 
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yfs = 500 GcV may border the 0(100 fb) level. This 
would translate into barely 5 x 10 4 events per 500 fb -1 
of integrated luminosity at the linac. While the former 
situation describes perhaps the most optimistic possibil- 
ity, cross sections of 0(10 fb) at y/s = 500 GeV should 
be quite common in a wide patch of the 2HDM parame- 
ter space, and correspondingly of 0(1 fb) at y/s = 1 TeV 
(see Table lTV]) , thus delivering rates that range from a few 
hundred to a few thousand events respectively. With this 
level of statistics (plus a comparable number of events 
from the H°A° channel) it should be perfectly possible to 
insure a comfortable tagging of the Higgs bosons in the 
clean environment of a linac machine. 

Most important, we have been able to demonstrate 
that the dominant quantum effects are, as we expected, 
primarily nurtured by the Higgs-mcdiatcd one-loop cor- 
rections to the A h° Z°/A° H° Z° vertices, and can be ul- 
timately traced back to the potentially enhanced 3H self- 
interactions Xsh (cf. Table II). This is indeed a trade- 
mark feature of the 2HDM, with no counterpart in the 
MSSM - in which SUSY invariance highly curtails the 
structure of the Higgs boson self-interactions and com- 
pels them to be purely gauge-like. Obviously, probing 
the structure of Higgs boson self-interactions - in the 
present case through the analysis of radiative corrections 
on direct production processes - is a most useful strat- 
egy to disentangle SUSY and non-SUSY Higgs physics 
scenarios. 

To be specific, for the analysis of the e + e~ —> 2H 
cross sections we have faced the computation of the full 
set of 0(a^, w ) quantum effects, among them the subset 
of 0{a 2 ew \\ H ) corrections, and supplemented them with 
the leading 0(a^ w ) pieces arising from the square of the 
one-loop diagrams involving Higgs boson self-couplings. 
These (finer) corrections are of 0(a 2 w X^ H ) (see, how- 
ever, below). Their inclusion is actually a must at large 
center-of-mass energies (i.e. once y/s has left behind the 
startup situation, and of course also in the 1 TeV regime). 
The reason is that the tree- level 0{a 2 ew ) and the one- 
loop 0(al w ) amplitudes virtually cancel out each other 
at high energies, owing to the characteristic large and 
negative quantum effects in this regime. We can con- 
vince ourselves of this fact from the severe depletion ex- 
hibited by the one-loop curves - see e.g. Figs. [T51ITT1 
- as compared to the tree-level ones. Typically, the ef- 
fect appears after we increase the center-of-mass energy 
by 100 GeV beyond the initial y/s = 500 GeV, and per- 
sists till 1 TeV and afterwards. At large energies, the 
inclusion of the 0{a 2 . w X\ H ) effects is a consistency re- 
quirement as it prevents the cross-section from becoming 
negative. These terms can actually be relevant also at 
the initial energy interval y/s = 500 — 600 GcV at which 
the linac will first operate; for example, for the scenarios 
with maximally enhanced 3H self-couplings, they may 
lead to an additional ~ 20% contribution to 6 r , hence 
pushing the overall quantum correction even higher. In 
this respect, we note that our earlier estimate of these 
higher order terms, Eq. (|52"]l . was able to hit the right 



order of magnitude, although it actually underestimates 
the maximum size revealed by the exact numerical anal- 
ysis - not too surprising from such a rough attempt at 
guessing their bulk size. 

Let us clarify that the two-loop diagrams involving 
the Higgs boson-mediated A° h° Z°/A° H° Z° vertex cor- 
rections, lead of course also to amplitudes of order 
0(a 2 w Xf H ). These are actually finite, and originate from 
the interference of the tree-level amplitude with the two- 
loop one. An obvious concern is then the following: can 
we safely neglect them? Upon inspection of the Higgs- 
mediated two-loop diagrams, and taking into account 
power counting and dynamical considerations, one can 
show that they are actually suppressed, and particularly 
so for those scenarios where cr^ 0+1 ^ (e + c~ — > 2H) is opti- 
mized. But, most significantly, there is a deeper obser- 
vation to be adduced at this point, which is well in the 
spirit of the effective field theory approach to QFT. The 
(finite) two-loop effects induced by the enhanced Higgs 
boson self-couplings can be conveniently reabsorbed in 
the form of one-loop corrections to the bare Higgs self- 
couplings, which first appear at the one-loop level in the 
process under consideration. We can iterate this algo- 
rithm at any order in perturbation theory by further re- 
absorbing the higher order Higgs mediated loops into the 
triple and quartic Higgs self-couplings introduced at one- 
loop order. In this way we can define a collection of effec- 
tive 3H and 4H self-couplings that encapsulate all these 
higher order effects. 

The above remark is an interesting one, as it tells us 
that the application of the stringent tree-level unitarity 
relations that we have used can be directly reiterated 
for these effective couplings. Therefore, no matter how 
big are the quantum corrections induced by the 3H self- 
couplings in subsequent orders of perturbation theory, 
their overall effect is constrained by the same set of for- 
mal unitarity relations. If expressed in terms of the orig- 
inal couplings, these relations will be equally stringent 
so that, on balance, the maximum enhancement capabil- 
ities of the 3H self- couplings are basically the same at 
any order of perturbation theory. This is the procedure 
adopted in the present work, and we believe it is per- 
fectly consistent, in the sense that we have retained the 
contributions to the e + e~ —> 2H amplitudes at leading 
order in the 3H self-couplings while employing the uni- 
tarity constraints of Ref. 67|-[6!|, also derived at leading 



order in X^h- 

How do our results compare with the expectations in 
the MSSM? As we have stressed in section II, a lot of work 
on Higgs boson production in e + e _ colliders has been re- 
ported in the literature thus far, but mainly within the 
context of the MSSM. In particular, supersymmetric ra- 
diative corrections to a(e + e _ — > h°A°) have been consid- 
ered in detail, although mostly in the context of LEP and 
including sometimes a more or less timid incursion into 
the TeV-class colliders, cf. Refs. [H [H [H [H M, HI ■ 
Let us try to see how some of these results compare with 
ours. For the sake of definitcness, let us concentrate on 
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the particular case of Ref. [37|, in which a(e + e~ — > A°h°) 
is analyzed in the MSSM for a linear collider operating at 
LEP 2 energy and further extended up to y/s = 500 GeV. 
In the most favorable regimes, the loop-corrected cross 
sections may reach up to barely 30 fb at y/s = 500 GeV, 
which falls in the ballpark of the results that we have 
obtained here within the 2HDM. This is not surprising 
because in both cases the leading amplitude (i.e. the 
lowest order one) is purely gauge and hence the order of 
magnitude of the cross-section is independent of whether 
we compute it in the MSSM or in a generic 2HDM. 

Quantum effects can be fairly efficient in the MSSM 
too, and they may also boost the tree-level predictions 
substantially (~ 20%). However, their size is never as 
big as the maximal 2HDM contributions; and, no less 
important for our conceptual understanding, the funda- 
mental origin of these quantum effects is completely dif- 
ferent from that of the 2HDM. As already mentioned 
in the introduction, the leading quantum effects in the 
MSSM are due to the genuine Higgs-quark, Higgs-squark 
and quark-squark-chargino/ncutralino Yukawa-like su- 
persymmetric couplings, whose contributions carry en- 
hancement factors of the fashion ~ nib tan /3/Myy, ~ 
m t (A t - fj,/ta,nP)/M% USY , m b (A b - n tan/3)/M% USY , 
none of them related to the structure of the Higgs poten- 
tial and hence none of them associated to Higgs boson 
self-couplings. 

In contradistinction to the MSSM, the 3H self- 
couplings in the 2HDM embody the full potential for 
triggering quantum effects in physical processes, and 
these effects become thence tied to the very structure 
of the Higgs potential. Even if quantitatively similar in 
some scenarios, quantum effects in SUSY and non-SUSY 
two-Higgs doublet extensions of the SM are, therefore, 
prompted by intrinsically different mechanisms. In prac- 
tice, the task of distinguishing whether the produced 
Higgs boson particles at a linac are supersymmctric or 
not can start by addressing the decaying signatures of the 
Higgs bosons into bb or r + r _ pairs, as sketched above. 
However, in certain regions of parameter space the num- 
ber of events can be by itself a direct signature of the sort 
of Higgs model that we have behind. For example, in 
the MSSM, the tree-level coupling C A o h o z o ~ cos(/3 — a) 
undergoes a severe suppression with growing values of 
tan j3 and/or Mao - as a depends on Mao in the MSSM. 
Thereby the predicted cross sections, either at the Born 
or at the one-loop level, are dramatically weakened as a 
function of Mao. For instance, cr' 0+1 '(e + e _ — > h°A°) at 
yfs = 500 GeV and M A o = 300 GeV may end up as tiny 
as - 0.01 fb (see e.g. Fig. 9 of Ref. [37|) while the 2HDM 
prediction is of O(10 fb), i.e. one thousand times bigger! 

But there is more to say here, as we have announced 
above. The identification strategy can be highly facili- 
tated through the interplay of further 2H (and eventually 
3H) channels, such as ([1]) and (pi), which are also highly 
distinctive in the 2HDM - see |27| and [1^]. In fact, for 
these channels the difference with the MSSM can be ap- 
parent already at the lowest order of perturbation theory. 



Let us be a bit more quantitative here. In Table [Xj we 
compare the numerical predictions on the cross-sections 
of these various processes in truly equitable conditions, 
i.e. for different sets of common parameters for all the 
processes, and for three realistic values of the center-of- 
mass energy of the ILC/CLIC (the latter operating al- 
ways at the highest edge of the energy band, although it 
could nominally reach the 3 TeV end) . 

To be more precise, the cross-sections that we are 
comparing on equal footing in Table [X] are the follo- 
wing: i) the one-loop corrected a^ +1 \e + e~ — > A°h°); ii) 
the leading-order 0{oq, w ) cross-sections of the two triple 
Higgs production processes e + e~ — > h H°A°/H n H + H~; 
and iii) cr(e+e" -> V*V* -> h°h° + X) at leading- 
order 0(a.g W ). Notice that the bulk of the contribution 
from processes ii) and iii) comes from the 0(a\ w \\ H ) 
and 0{a\ w \\ H ) parts, respectively. In all cases we take 
tan /3 = 1 , a = /3 and the maximum allowed value of | A5 1 
according to Table IVII1 that is, the most favorable sce- 
nario for h°A° pair production. We use Sets Till for the 
Higgs boson masses. The resulting predictions for the dif- 
ferent channels turn out to be highly illustrative of the 
complementary nature of such Higgs-boson production 
signatures. For a thorough discussion of the dynamical 
features and signatures of the multiparticle processes ii) 
and iii) , we refer the reader to the original references 
[27l [29j . Here we have recomputed the cross-sections ii) 
and iii) under the same kinematical conditions and in 
the very same region of the parameter space as the main 
processes i) under consideration, which is of course the 
indispensable first step that enables us to compare them 
all meaningfully. In this way, we are well prepared to 
infer the final strategy that can be designed to efficiently 
search for 2HDM Higgs bosons at the linear colliders. 

The strategy that follows from the results encapsulated 
in Table [X] and the information that we have gathered in 
the figures of the last section should be by now crystal- 
clear, and it comes down to the following two-step pro- 
cedure: 



• Step 1. At the vicinity of the startup energy 
of the ILC (approximately in the range s/s = 
500 ± 100 GeV), the exclusive neutral double Higgs 
boson channels e+e" -> A Q h°/A H - Eq. © - 
prove to be the dominant ones as compared to the 
triple Higgs boson production processes e + e~ — > 
3H and the inclusive double Higgs production ones 
e + e~ — >• 2H + X, see equations and ©. There- 
fore, in this lower energy band, the study of poten- 
tial signatures of new physics, and in particular the 
identification of the nature of the produced Higgs 
bosons at the ILC, must be conducted through the 
careful analysis of the quantum corrections affect- 
ing the two-body channels e + e~ — > 2H. Here is 
precisely where the detailed results of the present 
work could be most useful, and indeed should be 
the first analysis to be implemented when the lin- 
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Process 


|cr(^ = 0.5TeV)fb 


a 


(l.OTeV)fb 


<r(1.5TeV)fb 


Set I 


h u A u 


34.13 






2.89 


0.70 


h°H°A° 


3.09 






8.58 


5.17 


H°H+H _ 


6.75 






19.65 


12.27 


h°h° + X 


13.29 






79.00 


146.08 


Set II 


h u A u 


26.71 






4.07 


1.27 


h°H°A° 


0.02 






5.03 


3.55 


H°H+H- 


0.17 






11.93 


8.39 


h°h° + X 


1.47 






17.36 


38.01 


Set III 


h u A u 


11.63 






6.11 


2.52 


h°H°A° 


below threshold 




1.25 


1.33 


H°H + H~ 


below threshold 




0.69 


2.14 


h°h° + A 


0.92 






9.72 


23.40 



TABLE X: Comparison of the predictions for the cross sec- 
tions (in fb) corresponding to some of the Higgs boson pro- 
duction processes fTJ, ([2j and (|3|. The results are obtained 
for the Sets I, II and III of Higgs boson masses (cf. Ta- 
ble [IVJ) tan/3 = 1, a = j3, and three different values of the 
center-of-mass energy. We observe a great complementarity 
between the different channels at different energies: the exclu- 
sive 2H channels <(3j are dominant at the ILC startup energy 
y/s — 0.5 TeV whereas the others dominate at higher ener- 
gies, specially the inclusive ZR + X one {TJ- which is triggered 
mainly by weak gauge boson fusion. 

car colliders arc set to work in the future; 

• Step 2. At higher energies, however, say for 
y/s > 600 GeV and, for that matter, in the en- 
tire TeV range (hence at the highest nominal en- 
ergy regimes scheduled for the ILC, and specially 
for CLIC), the influence of the Z -propagator sti- 
fles dramatically the cross-section of the exclusive 
neutral double Higgs boson channels e + e~ — > 2H. 
Therefore, in these high energy domains, they can 
no longer compete with the mechanisms providing 
multiparticle final states, such as ([l} and ([2]). The 
latter process, e + e~ — > V*V* — > hh + X, becomes 
indeed the most efficient one at the highest energies 
since it is not crippled by the s-channel propaga- 
tor; quite on the contrary, its cross-section increases 
steadily with the energy (see Ref. [29[ for details). 
Therefore, in the upgraded phase of the ILC, and 
certainly for the CLIC collider, one must concen- 
trate all the search strategy power on looking for an 
anomalously large number of inclusive Higgs boson 
pairs of the type h°h°, H°H°, h°A° and A A , which 
should emerge mostly acollinear and dynamically 
thrustcd along the beam direction, rather than ap- 
pearing in a simple back-to-back geometry charac- 
teristic of the two-body process e + c~ — > 2H. Notice 
furthermore that, in contradistinction to the latter, 
the fused pairs consist of identical Higgs bosons. 
Unmistakcably, the predicted numbers in Table [X] 
could not, by any means, be accounted for if they 



were to be ascribed to a supersymmetric origin. Fi- 
nally, the cross-correlation of these higher energy 
effects with the previous ones at lower energies, 
which we could track very well while upgrading the 
ILC collider from y/s = 0.5 TeV up to y/s = 1.5 TeV 
- and eventually till yfs = 3 TeV (through CLIC) - 
should provide plenty of unambiguous evidence of 
non-supcrsymmctric Higgs physics 14 . 

Let us finally introduce a remark that hints at the po- 
tentially far reaching implications of the 2HDM dynamics 
in different sectors of Particle Physics. The triple Higgs 
boson self-couplings could also have an indirect impact on 
the best high precision observables at our disposal, even 
well before the ILC/CLIC colliders are commissioned and 
put effectively to work. For example, they could have a 
bearing on the famous electroweak precision parameter 
Ar - see Eq. (|llj) . Since our main aim here has been to 
exploit the leading contributions from the Higgs boson 
self-couplings in the arena of the direct production pro- 
cesses ([3]), the Ar part did not enter our quantum com- 
putation; it would only enter at two-loop level. However, 
we suspect that the influence of the 3H self-couplings 
on Ar might have a real interest per se, for it could in- 
duce a correction to 5p and ultimately trigger a shift in 
the W mass, and the result might well compete with the 
highly accurate calculations that have been performed on 
this observable within the alternative framework of the 
MSSMpjj]. 

To summarize, if the genuine enhancement proper- 
ties of the 2HDM effectively hold in the real world, the 
combined analysis of the exclusive double Higgs produc- 
tion (2H) processes at the startup energy of the ILC, 
= 500 GeV, and subsequently of the triple (3H) and 
inclusive double Higgs production processes (2H + X) 
at the upgraded ILC/CLIC (y/s > ITcV), might reveal 
strong hints of Higgs boson physics beyond the SM. The 
bottom line of our study is that the physics of the linear 
colliders can be truly instrumental to unveil the nature 
of the Higgs bosons. These bosons may, or may not, have 
been discovered during the LHC era while the linear col- 
liders take the floor, but the real issue at stake here is 
of highest priority and should strengthen the need for 
such machines. High precision Higgs boson physics in a 
linac can indeed provide a keen insight into the most sen- 
sitive building blocks of modern gauge theories of weak 
and electromagnetic interactions, viz. in the very core 
architecture of the Electroweak Symmetry Breaking. 



Another production mechanism where the outstanding enhan- 
cing properties of the general 2HDM Higgs-boson self-couplings 
can be tested - and correlated with the processes under consi- 
deration - is the one-loop production of neutral Higgs bosons in 
association with the Z gauge boson: e + e — — > Z°h (h = h° , H°), 
also called "Higgs-strahlung" processes; see [89l ] for a thorough 
study at the quantum level in the general 2HDM. 
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